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HEAT  TRANSFER  BY  FREE  CONVECTION  FROM  HORIZONTAL 
CYLINDERS  IN  DIATOMIC  GASES* 
By  R.  Hermann 

I.   DETERMINATION  OF  HEAT -TRANSFER  LAW  FOR  HORIZONTAL  CYLINDER  FROM 

TESTS  WITH  PARTICULAR  ACCOUNT  TAKEN  OF  THE 

TEMPERATURE  CHARACTERISTIC  Te^ 

1.  Introductory  Remarks 

The  case  of  the  horizontal  cylinder  is  of  particular  importance  in 
the  study  of  heat  transfer  by  free  convection  for  the  following  reasons: 
in  the  first  place,  next  to  the  rectangular  plate  it  represents  the 
simplest  two -dimensional  case;  and  second,  a  very  wide  range  of  measure- 
ments is  possible,  from  the  finest  electrically  heated  glow  lamp  wires 
to  pipes  heated  by  liquids  or  gases  flowing  through  them. 

To  investigate  free-convection  flow  from  the  point  of  view  of 
similarity  considerations,  it  is  convenient  to  consider  the  case  of 
small  temperature  differences  between  the  heated  body  and  the  surround- 
ings; in  this  case  all  the  properties  of  the  medium,  even  the  density 
(ref.'l,  p.  429),  in  the  entire  temperature  field  may  be  assumed 
constant.   The  case  of  large  temperature  differences,  for  which  the 
variation  in  properties  over  the  temperature  range  can  no  longer  be 


"Warmeiibergang  bei  freier  Stromung  am  wagrechten  Zyllnder  in 
zweiatomigen  Gasen."  VDI  Forschungsheft,  No.  379,  1936,  pp.  1-24. 

-'-Accepted  as  dissertation  by  the  Technical  High  School  at  Aachen 
with  the  approval  of  Prof.  Dr.-Ing.  C.  Wieselsberger  and  Prof.  Dr. 
W.  Muller.   The  tests  were  carried  out  in  the  division  for  applied 
mechanics  and  thermodynamics  of  the  Physical  Institute  at  the  University 
of  Leipzig.   The  author  takes  this  occasion  to  thank  Prof.  Dr.  L. 
Schiller  for  the  interest  which  he  took  in  this  work  and  for  his  valua- 
ble advice. 
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neglected,  is  treated  separately.   As  first  pointed  out  by  W.  Nusselt 
(ref.  2),    it  follows  from  the  differential  equations  of  free  convection, 
together  with  the  uniquely  determining  boundary  values  (ref.  1,  p.  428) 
for  the  case  of  small  temperature  differences,  that  the  nondimensional 
heat-transfer  parameter  Nu  (Nusselt  number^)  is  determined  by  the  lift 
coefficient  Gr  (Grashof  number)  and  the  characteristic  Pr  (Prandtl 
number)  for  the  molecular  constitution  of  the  gas'^: 

Nu  =  F(Gr,  Pr)  (l) 

where 

Nu  =  od/X,     Gr  =  d3  g  (3  0/v  2   Pr  =  v/a  (2) 

a  is  the  heat-transfer  coefficient  (cal/cm^) (sec) (  C),  d  is  the  cylin- 
der diameter  (cm),  g  is  acceleration  of  gravity  (cm/sec^),  X  is  the 
heat  conductivity  (cal/(cm) (sec) (°C) ) ,  v   is  the  kinematic  viscosity 
(cm^/sec),  a  is  the  temperature  conductivity  (cm^/sec),  p  is  the  coef- 
ficient of  expansion  (°C~  ),  and  ®  is  the  temperature  difference  (°C) 
between  body  (t^)  and  medium  (t„) . 

The  case  of  large  temperature  difference  between  the  body  and  the 
surroundings,  for  which  the  gas  properties  vary  over  the  field,  was 
similarly  first  considered  by  W.  Nusselt  (ref.  2).   The  following  ex- 
pression is  obtained  for  the  heat  transfer  under  the  restricting  con- 
dition that  the  temperature  dependence  on  v   and  a  may  be  repre- 
sented by  exponential  laws  in  the  absolute  temperature,  the  exponents 
of  which  for  gases  of  the  same  substance  must  be  equal'^; 

Nu  =  F(Gr,  Pr,  Te)  (3) 

The  form  here  chosen  for  the  temperature  characteristic^ 


■^In  the  notation  of  the  characteristic  numbers,  the  proposals  of 
the  Heat  Conference  at  Koln  in  1931  are  followed.   See  Forschg.  Ing.-Wes., 
Bd.  2,  p.  380,  1931. 

^See  eq.  (5). 

■^This  is  approximately  the  case,  for  example,  for  air,  oxygen,  and 
hydrogen;  see  reference  1,  p.  433. 

^The  abbreviation  Te  is  chosen  to  conform  to  the  other  abbrevia- 
tions for  the  characteristic  numbers  and  to  remind  the  reader  that  it 
refers  to  a  temperature  relation. 
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Te  =  0/T„  (4) 

(T  ,  Too^  (°K)  =  absolute  temperatures  of  the  body  and  the  surroundings, 

respectively)  differs  inessentially  from  the  nondimensional  ratio 
T^/T«,  chosen  by  Nusselt,  namely  by  the  constant  -1.   Equality  of  Gr, 

Pr,  and  Te  in  two  different  cases  means  complete  similarity,  that  is, 
similarity  of  the  velocity  field,  the  temperature  field,  and  the  field 
of  all  the  coefficients  characteristic  of  the  substance  (ref.  1,  p.  434), 
and  therefore  also  equality  of  the  Nu  number.   The  possibility  still 
remains,  however,  that  similarity  may  exist  at  the  same  time  that  the 
previously  mentioned  condition  of  the  power  law  is  not  satisfied. 

In  the  present  report  the  characteristics  represented  in  equa- 
tion (3)  will  be  computed  for  the  heat  transfer  from  heated  wires  and 
pipes  in  air,  hydrogen,  and  oxygen;  and  from  these  values  the  most 
probable  form  of  the  heat-transfer  law,  obtainable  at  present,  will  be 
determined  for  the  case  of  the  horizontal  cylinder  in  diatomic  gases 
(Pr  =  0.74).^^  The  dependence  of  Nu  on  the  temperature  characteris- 
tic Te,  which  follows  from  the  similarity  theory  and  is  illustrated  by 
an  example  (ref.  1,  fig.  l)  as  a  third  independent  variable,  will  be  ver- 
ified for  several  tests  in  the  range  of  Gr  from  10"'^  to  10.  For  the 
smallest  values  of  Gr,  the  effect  of  Te  considerably  outweighs  that 
of  Gr.   In  the  region  of  large  Gr  from  10'^  to  10',  on  the  contrary, 
Te  is  practically  without  effect.   The  large  scatter  of  the  test  points 
in  the  nondimensional  representation  of  the  test  results  of  Wusselt 
(ref.  2)  and  Davis  (ref.  3)  (±32  percent  in  the  case  of  Nusselt)  is  due 
to  the  fact  that  the  parameter  Te  is  not  taken  into  account  in  accord- 
ance with  equation  (l).   The  dependence  of  Nu  on  Gr,  Pr,  and  Te 
will  be  theoretically  clarified  hereinafter.   The  effect  of  Te  for 
small  Gr  may  also  be  correctly  estimated  quantitatively. 

From  the  tests  conducted  by  various  investigators  for  different 
finite  Te  values,  the  limiting  law  of  small  temperature  difference  is 
determined  in  each  case  by  extrapolation  to  Te  =  0.   This  is  of 
primary  theoretical  and  practical  significance  since  it  is  identical 
with  the  theoretically  simpler  case  (l)  and  is  thus  free  from  the  pre- 
viously mentioned  restrictive  assumption  with  regard  to  the  temperature 
dependence  of  the  constants  defining  the  properties  of  the  substance. 


For  the  values  of  X   and  v   required  for  computing  the  char- 
acteristics, values  are  assumed  for  the  temperature  t^  (in  agreement 

with  E.  Schmidt).   The  similarity  consideration  leaves  the  choice  of 
the  reference  temperature  free. 
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2.  Radiation 

From  the  measured  over-all  heat  losses  directly  determined  in  the 
testSj  the  amount  due  to  radiation  must  be  subtracted.   For  this  pur- 
pose, a  well  known  eqtiation  (ref .  4,  eq.  38,  p,  232)  was  employed;  the 
best  mean  values  at  present  known  were  substituted  as  radiation  con- 
stants for  the  metal  cylinder  surfaces,  as  follows:   for  platinum,  nickel, 
silver,  and  tantalum,  the  values  given  by  H.  Schmidt  (ref.  5);  and  for 
copper  and  iron,  the  values  given  by  Grober  (ref.  6,  p.  196).   Small 
deviations  as  compared  with  the  values  used  by  the  investigators  them- 
selves are  without  significance  in  the  range  of  small  values  of  Gr 
because  the  radiation  component  in  the  case  of  thin  wires  amounts  to 
only  a  few  percent  of  the  heat  of  convection.   For  the  tests  of  Wamsler 
and  of  Koch  on  thick  pipes,  however,  for  which  the  radiation  losses  are 
up  to  60  percent,  an  accurate  knowledge  of  the  radiation  coefficients 
is  required.   These  coefficients  were  obtained  by  Wamsler  (ref.  13) 
through  his  own  radiation  tests.  After  correct lon^  of  computations  of 
reference  13,  the  values  obtained  herein  for  the  radiation  heat  are 
0.2  to  2  percent  larger. 


3.  Property -Determining  Constants  of  Gases 

In  order  to  compute  Nu,  Gr,  and  Pr,  values  of  the  following 
properties  are  required:  density  p,  dynamic  viscosity  n,  heat  con- 
ductivity X,  and  specific  heat  c  for  air,  hydrogen,  and  oxygen. 

(a)  Constants  for  air.  -  The  density  was  computed  throughout  the 
temperature  and  pressure  range  according  to  the  ideal  gas  law^ . 

The  true  specific  heat  for  constant  pressure  was  obtained  from  an 
equation  given  by  Holborn  and  JakobS  for  the  mean  specific  heat. 

The  temperature  dependence  of  the  dynamic  viscosity  was  determined 
by  graphical  adjustment  of  the  values  given  by  Erk^.   These  values 


^Wamsler  (ref.  13)  computes  incorrectly  with  the  previously  men- 
tioned radiation  equation  without  the  factor  of  the  area  ratio  in  the 


denominator . 

7 

'Landolt-Bornstein:  5th  ed. ,  vol.  I,  p.  43. 

°Landolt-Bornstein:  5th  ed.,  vol.  II,  p.  1274. 

^Landolt-Bornstein:  5th  ed.,  vol.  I,  pp.  177-181,  Eg.  I,  pp.  143- 
144,  Eg.  lla,  pp.  138-141. 
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differ  only  slightly  (at  most,  4  percent)  from  the  graphical  mean 
values  given  by  Erk  from  0°  to  700°  C-^"-". 

The  heat  conductivity  at  high  temperatures  was  obtained  by  a  gas 
kinetic  extrapolation,  since  the  test  values  collected  by  Jakob   lie 
only  between  0°  and  212°  C.   The  Nusselt  formula  (ref.  1,  p.  491; 
ref.  6,  p.  192),  which  deviates  from  the  available  test  values  by  at 
most  2.8  percent,  is  employed  herein. 

From  the  gas  kinetic  relation 

Pr  =  Id  c/X  =  (n  +  2)/(n  +  4.5)  (5) 

where  n  is  the  number  of  degrees  of  freedom  of  the  molecular  motion 
according  to  which  the  Pr  number  of  a  gas  is  independent  of  the  tem- 
perature (and  moreover  is  the  same  for  all  gases  with  the  same  number 
of  atoms);  there  is  obtained  for  air  a  small  decrease  of  Pr  =  0.739 
for  0°  C  to  0.721  at  200°  C  and  0.679  at  1000°  C  (value  obtained  from 
formula  for  n=5  is  Pr  =0.737),  which,  however,  is  of  no  signifi- 
cance because  it  lies  within  the  uncertainty  limits  of  the  jj.  and   X 
values.   The  same  holds  for  the  differences  in  the  Pr  numbers  among 
air,  hydrogen,  and  oxygen. 

(b)  Constants  for  hydrogen  and  oxygen.  -  For  the  kinematic  vis- 
cosity and  heat  conductivity  the  values  given  by  Davis  (ref.  3)  are 
employed,  which  also  for  air  are  in  good  agreement  (deviation,  at  most 
2.7  percent)  with  the  chosen  values.   A  knowledge  of  the  specific  heat 
is  unnecessary  if  the  Pr  value  for  0°  C  is  satisfactory;  the  values 
are  0.717  for  hydrogen  and  0.731  for  oxygen-^^. 

(c)  Gas  constants  for  high  pressures.  -  From  the  tests  of  Petavel, 
which  extend  up  to  160  atmospheres,  only  those  up  to  40  atmospheres 
were  computed  in  order  that  reliable  gas  constants  would  be  available. 
In  this  pressure  range  and  up  to  1000°  C,  the  computation  was  conducted 


-'-'\ien-Harms :  Handb.  d.  Exp.  Phys.,  IV,  4,  p.  531. 

Landolt-Bbrnstein:  5th  ed.,  vol.  II,  p.  1304. 

A.  Busemann  in  Wien-Harms:  Handb.  d.  Exp,  Phys.,  IV,  1,  p.  359. 
This  relation  may  be  obtained  from  an  equation  given  by  Eucken  (Physik 
Z.,  vol.  14  (1913),  p.  324. 

"I  r! 

A.  Busemann  in  Wien-Hanns:  Handb.  d.  Exp.  Phys.,  IV,  1,  p.  362. 
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with  the  normal  values  since  the  dependence  on  the  pressure  of  X  and 
\x     and  the  deviations  of  p  from  the  ideal  gas  law  are  either  not  known 
at  all  or  are  only  partly  known,  and  the  known  deviations  amount  to  only 
a  few  percent-'-^. 


4.  Individual  Investigations-'-^ 

(a)  General.  -  The  nondimensional  representation  of  the  results  is 
effected  according  to  equation  (3)  hy  logarithmic  graphs  with  Nu  as 

a  function  of  Gr  and  with  Te  as  the  only  parameter,  since  Pr 
(according  to  sec.  3a)  may  he  considered  as  constant.   Even  if  the 
slight  variation  of  Pr  there  indicated  is  considered  as  real,  there 
would,  according  to  a  test  result  of  Davis  (ref .  3)  on  the  effect  of 
the  Pr  number,  he  obtained  only  displacements  of  the  curves  of  the 
order  of  magnitude  of  l/lOO  of  those  which  can  actually  be  observed. 

In  table  I  are  summarized  the  test  objects  and  test  conditions  of 
the  individual  investigators  and  the  range  of  Te  used.   In  order  to 
judge  to  what  extent  the  condition  of  infinite  extension  of  the  fluid 
was  satisfied,  the  ratio  of  the  height  H  of  the  surrounding  space  to 
the  cylinder  diameter  d  is  shown.   In  figures  1  to  7,  approximately 
equal  Te  values  are  indicated  by  the  same  symbols.   The  decrease  in 
Nu  with  increasing  Te  is  particularly  evident  in  the  individual 
point  groups  (because  of  the  same  wire),  whereas  the  values  for  different 
wires  naturally  show  larger  scattering  as  compared  with  one  another. 
The  curves  Te  =  constant  were  obtained  by  a  graphical  adjustment  pro- 
cess and,  similarly,  the  extrapolation  points  and  curves  for  Te  =  0 
and  the  dependence  of  Ru  on  Te  (relative  decrease  in  Nu  for 
ATe  =  l).   The  indicated  minimum  and  maximum  values  hold  for  the  indi- 
vidual wires  (point  groups);  the  indicated  mean  values,  for  the  individ- 
ual investigators. 

(b)  Results  of  Ayrton  and  Kilgour  (ref.  l),  figure  1.  -  The  authors 
of  reference  1  give  in  the  form  of  diagrams  for  each  wire  the  total  heat- 
transfer  coefficient  as  a  function  of  the  wire  temperature;  two  wires 
(0.0206  and  0.0282  cm  diam. )  because  of  the  very  large  scattering  and  a 
third  wire  (0.0102  cm  diam.)  because  of  its  strong  deviation  in  position 
and  inclination  were  excluded  from  the  evaluation  of  the  tests.   The 


^^Landolt-Bornstein:  5th  ed. ,  Eg.  I,  p.  64  ff,  Eg.  Ila,  p.  144. 

Acknowledgement  is  made  to  Dr.  K.  Winkler,  Leipzig,  for  his 
assistance  with  the  computations. 
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points  computed  from  the  low  temperatures  of  40°  to  80°  C  were  simi- 
larly not  used  on  account  of  the  large  scattering. 

(c)  Results  of  Langmuir  (ref.  8),  figure  2.  -  Langmuir  gives  in  a 
tatle  the  adjusted  values  of  the  total  energy  losses  as  a  function  of 
the  absolute  wire  temperature  increasing  in  stages  of  200°  C.   Only 
the  tests  up  to  1300°  K  were  computed  herein  in  order  to  have  reliable 
values  available  for  v   and   X.  The  convergence  of  the  curves 
Te  =  1.33  and  Te  =  3.33  shows  that  the  dependence  of  Nu  on  Te 
decreases  with  increasing  Gr,  the  decrease  independently  determined 
for  the  individual  wires  being: 


log  Gr 

Decrease,  percent 


0.1-4 
7.2 


0.8-4 
7.5 


0.6-3 
5.4 


0.5-2 
3.2 


0.4-1 
2.2 


This  experimental  confirmation  of  the  results  following  from  theore- 
tical considerations  (sec.  6b)  acquires  significance  in  that  the  tests 
of  Langmuir  are  characterized  by  a  large  change  of  the  Te  number  and 
smallest  scattering.   The  extrapolation  from  Te  =  1  to  Te  =  0,  on 
account  of  the  absence  of  test  points  for  this  value,  was  effected  with 
the  mean  value  of  13.2  percent  taken  from  reference  7. 

Langmuir  (ref.  9)  also  discusses  his  own  tests  at  various  air 
pressures  (lO  to  760  mm  Hg)  and  variable  room  temperature  (t<»  =  -190°  C 
to  600°  C)  which  on  account  of  the  strong  change  in  the  expansion  coef- 
ficient (p  =  1/Too)  would  be  of  great  interest  and  the  only  tests  of 
their  kind.   Unfortunately,  these  tests  have  not  as  yet  been  published, 
and  the  test  data  are  not  available  in  a  form  in  which  they  may  be 
evaluatedlS. 

(d)  Results  of  Bijlevelt  (ref.  lO),  figure  3.  -  Since  the  tests  of 
Bijlevelt  were  primarily  for  the  purpose  of  investigating  forced  con- 
vection, the  experimental  setup  was  designed  for  the  sensitivity  re- 
quirements of  this  flow.   This  explains  the  greater  scatter  of  the  test 
points  of  reference  10  as  compared  with  those  of  other  investigators. 
For  each  individual  wire  (with  the  exception  of  the  Ta  and  Ni  wire), 
a  grouping  is  nevertheless  observed  with  respect  to  Te  in  the  usual 
sense. 

(e)  Results  of  Kennelly,  Wright,  and  Bijlevelt  (ref.  ll),  fig- 
ure 4.  -  The  investigation  of  reference  11  is  noteworthy  in  that  through 
the  pressure  (p, v)  changes,  even  for  constant  wire  diameter  and  constant 
surrounding  temperature,  a  change  is  effected  in  the  Gr  number  so  that 


-1  (2 

For  these  data  the  author  is  indebted  to  a  personal  communication 
from  Dr.  I.  Langmuir. 
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curves  of  Te  =  constant  are  here  obtained  from  tests  on  a  single 
wire.   The  nondimensional  representation  of  these  tests  (fig.  4)  has 
been  given  previously  (ref.  l) .   The  authors  give  for  each  wire  a  dia- 
gram with  the  heating  current  as  a  function  of  the  pressure  for  three 
or  four  different  temperature  differences,  as  parameter  of  which  the 
values  for  five  different  pressures  were  computed  in  terms  of  the 
characteristics  and  were  connected  by  the  curves  Te  =  constant.   The 
extrapolation  to  Te  =  0  shows  extremely  good  agreement  for  both  wires. 
The  values  of  the  wire  of  medium  diameter  (d  =  0.02616  cm)  were  not 
plotted  because  of  the  irregularities  obtained. 

(f)  Results  of  Petavel  (ref.  12),  figure  5.  -  From  the  over-all 
heat-transfer  coefficients  given  for  constant  pressure  as  a  function 
of  the  temperature  difference,  values  were  computed  for  200°,  400°, 
600°,  800°,  and  1000°  C  temperature  differences  up  to  pressures  of 

40  atmospheres.   For  the  sake  of  clarity,  however,  only  the  curves  for 
the  extreme  Te  values  are  plotted  in  figure  5.   Since  the  heat  con- 
ductivity and  density  of  hydrogen  differ  from  the  values  of  air  by  the 
factors  6  and  14,  respectively,  these  tests  would  have  been  particu- 
larly suitable  for  checking  the  similarity  law.   Unfortunately,  however, 
these  tests  are  evidently  unreliable  because  of  the  far  too  small 
jacket  pipe  (H/d  =  18.6  as  the  smallest  value  of  all  investigators). 
The  similarity  law  for  the  three  gases  is  confirmed  only  for  Te  =  0.69 
and  Gr  >  10^.   In  the  case  of  air  and  oxygen,  there  is  a  splitting 
with  respect  to  the  dependence  on  Te,  but  in  the  reverse  sense  from 
that  observed  for  all  wires  of  the  other  investigators.   In  the  case 
of  hydrogen,  a  dependence  on  Te  in  the  correct  sense  is  observable 
only  at  the  smallest  values  of  Gr.   Finally,  the  Nu  values  also  for 
small  Gr  lie  too  high  as  compared  with  the  values  of  all  other  in- 
vestigators and  the  rise  in  the  Nu  with  increasing  Gr  has  too  great 
a  lag.   For  these  reasons,  the  results  of  Petavel  for  determining  the 
heat-transfer  law  were  not  used. 

(g)  Results  of  Wamsler  (ref.  15),  figure  6.  -  Wamsler's  curve, 
which  was  obtained  on  the  basis  of  correctly  computed  radiation  loss 
(sec.  2),  shows  a  dependence  on  Te  such  that  the  mean  values  for 

Te  =  0.60  lie  about  5  percent  higher  than  those  for  Te  =  0.20.   This 
small  splitting  of  the  effect  of  Te  which  occurs  in  the  opposite 
sense  to  that  in  the  case  of  small  Gr,  however,  in  all  probability 
does  not  contradict  the  facts  but  may  be  explained  by  the  consideration 
that  the  radiation  coefficients  in  the  entire  temperature  range  from 
50°  to  270°  C  were  assumed  constant,  whereas  in  general  for  metals  there 
has  been  established  an  increase  of  the  radiation  coefficient  with  the 
temperature.   If  in  the  previously  mentioned  temperature  range  a  rise 
in  the  radiation  coefficients  of  about  5  percent  is  assumed,  as  follows 
for  cast  iron  from  Wamsler's  determination  of  the  radiation  coefficient 
and  for  wrought  iron  from  the  determination  of  Wusselt  (ref.  14),  the 
observed  splitting  of  the  effect  of  Te  is  neutralized  and  the  curve 
shown,  valid  for  all  Te  values,  is  obtained.   The  test  values  of  the 
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copper  pipe  and  of  a  wrought  iron  pipe  (of  5.9-cm  diam. )  were  not  em- 
ployed in  figure  6  'because  of  the  obvious  irregularities  indicated  by 
the  scattering. 

(h)  Results  of  Koch  (ref.  15),  figure  7.  -  The  nondimensional 
characteristics  were  redetermined  herein  with  the  values  X   and   v 
for  wall  temperature,  since  the  characteristics  given  by  Koch  were  com- 
puted for  a  mean  temperature  of  the  substances.   The  test  results  show 
an  unusually  small  scattering  which  is  to  be  ascribed  a  careful  account- 
ing for  various  factors  (e.g.,  determination  of  the  temperature  of  the 
air  in  the  room  and  at  the  walls,  and  also  the  end  effect  of  the  pipes). 
For  the  two  intermediate  pipes,  after  graphical  adjustment,  a  decrease 
of  Nu  for  Te  =  0.07  to  0.52  by  about  2.5  percent  is  observed,  where- 
as in  the  case  of  the  smallest  pipe  no  ordering  is  observed  and  in  the 
case  of  the  largest  pipe  there  is  a  partial  decrease  and  a  partial 
scattering.   In  agreement  with  the  results  of  Wamsler,  the  dependence 
on  Te  is  thus  practically  zero  also  for  large  Gr   (lO  to  10  ) . 


5.  Summary  of  Test  Results 

The  results  of  the  different  investigators  as  regards  their  agree- 
ment may  be  compared  in  two  groups.   In  the  region  of  small  Gr 
(10""^  to  10),  the  tests  of  Ayrton  and  Kilgour,  Langmuir,  Bijlevelt,  and 
Kennelly  and  coworkers  are  in  agreement.   From  these  investigations  are 
determined,  on  the  one  hand,  the  curves  Te  =  0  for  the  limiting  law 
of  small  temperature  difference  (fig.  8)  and,  on  the  other  hand,  the 
curves  Te  =  0.65  corresponding  to  an  intermediate  test  range  (fig.  9). 
For  Te  =  0  the  maximum  scatter  of  the  results  of  the  four- investigators 
is  10.5  percent,  the  minimum  5  percent,  the  corresponding  values  for 
Te  =  0.65  being  22  and  8.5  percent,  respectively.   The  middle  curves 
of  figures  8  and  9  obtained  for  Te  =  0  and  Te  =  0.65  give  for  this 
change  in  Te  an  average  change  in  Nu  of  14.7  percent  (fig.  10). 
A  relative  decrease  in  the  Nu  number  by  22  percent  with  increase  of 
Te  from  0  to  1  is  thus  obtained  as  the  experimental  mean  value  of  the 
effect  of  Te  for  small  Gr. 

In  the  region  of  large  Gr  (lO  to  10'),  the  test  values  of 
Wamsler  and  Koch  from  energy  measurements  and  the  values  obtained  from 
the  temperature  field  measurements  of  Jodlbauer  (ref.  16)  by  integrat- 
ing over  the  cylinder  perimeter  are  shown  in  figure  10.   As  follows 
from  the  discussion  of  the  results  of  Wamsler  and  Koch,  the  values  of 
Koch  must  be  considered  as  more  reliable.   There  is  nevertheless  an 
uncertainty  in  the  values  due  to  the  uncertainty  in  the  radiation  com- 
ponents which  lie  between  40  and  60  percent.   On  the  other  hand,  the 
temperature  field  measurements  of  Jodlbauer  are  free  from  this  source 
of  error.   Since,  moreover,  the  results  of  reference  16  are  in  good 
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agreement  (mean  deviation,  4  percent)  with  the  theoretical  solution 
given  in  Part  II  they  were  assigned  a  greater  weight  in  the  evaluation. 
Jodlhauer  himself,  however,  considers  his  energy  measurements,  which 
are  about  10  percent  higher  and  lie  near  the  values  of  Koch,  to  be 
correct  and  his  field  measurements  in  error  because  of  the  introduction 
of  the  thermocouple.   The  test  points  of  Wamsler  are,  on  the  average, 
18  percent  higher  than  the  theoretical  values,  and  those  of  Koch, 
11  percent  higher.   The  interpolation  curve  is  not  everywhere  satis- 
factory as  regards  its  slope  and  curvature  and  could  be  determined  only 
with  a  partial  deviation  from  the  test  values. 

The  differences  still  existing,  after  extensive  adjustment,  among 
the  individual  investigators  and  the  difficulty  of  interpolation  show 
that  the  heat-transfer  curve  is  still  not  conclusively  determined  with 
the  desirable  accuracy  of  several  percent,  but  that,  on  the  contrary, 
further  tests  are  required.   The  given  interpolation  curve  nevertheless 
represents  the  most  probable  curve  of  the  heat -transfer  law,  obtainable 
at  the  present  time,  for  diatomic  gases  with  account  taken  of  the  de- 
pendence on  Te  in  the  range  of  Gr  from  10"  to  10'.   The  niimerical 
values  of  this  curve  are  given  in  table  II. 


6.  Qualitative  Theoretical  Interpretation  of  Heat-Transfer  Law 

(a)  Streamlines.  -  It  is  necessary  to  explain  theoretically  the 
dependence  of  the  Nu  number  on  Gr  and  Te,  determined  experimentally 
in  sections  4  and  5  (fig.  10 ),  and  also  the  dependence  on  Pr  as  ob- 
tained by  A.  H.  Davis  (ref .  3)  for  liquids.   A  solution  of  the  differ- 
ential equations  under  the  assumption  of  a  thin  (as  compared  with  the 
cylinder  diam. )  heated  laminar  layer  is  given  in  Part  II  for  diatomic 
gases  (Pr  =  0.74)  and  small  temperature  differences  (i.e.,  approximately 
constant  properties,   Te  «  O)  and  is  valid  for  mean  values  of  Gr  of 
about  10"^  to  3>:i0®.   On  the  basis  of  these  theoretical  results  and  with 
the  aid  of  suitable  assumptions,  the  dependence  of  Nu  on  small  and 
large  Gr  and  also  on  Pr  and  Te  will  be  discussed  theoretically 
without  any  necessary  implications  of  finality  in  the  conclusions. 

The  streamline  plot  (fig.  11)  obtained  from  the  solution  and  the 
isotherms  of  figure  12  show  that  cold  air  streams  from  the  bottom  and 
sides  of  the  cylinder  with  increasing  velocity  in  laminar  flow  upwards 
along  its  surface,  is  deflected  away  from  the  surface  in  the  region  of 
the  upper  stagnation  point,  and  forms  a  rising  current  of  warm  air. 
With  increasing  Gr  of  the  cylinder,  the  thickness  of  the  streaming 
layer  relative  to  the  cylinder  diameter  decreases.   Since  the  velocities 
increase  at  a  higher  rate,  however,  the  Reynolds  number  Re  of  the 
boundary  layer  increases  so  that  the  boundary  layer  finally  becomes 
turbulent.   According  to  the  schlieren  photographs  of  figure  23,  this 
occurs  in  the  neighborhood  of  the  upper  stagnation  point  at  Gr  =  3.5x10  . 
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With  further  Increasing  Gr,  the  point  of  turbulence  transition  travels 
upstream  (fig.  26)  and  reaches  the  equatorial  region  at  Gr  =  3^10  . 
With  decreasing  Gr  of  the  cylinder,  the  Re  of  the  "boundary  layer 
likewise  decreases,  and  at  the  same  time  the  thickness  of  the  layer  in 
relation  to  the  cylinder  diameter  increases  so  that  the  given  solution 
of  the  differential  equations  finally  (below  Gr  =  10^)  becomes  in- 
valid.  At  very  small  Gr  (very  small  Re) ,  the  isothermal  picture  of 
the  free  flow  in  the  neighborhood  of  the  body  becomes  increasingly 
similar  to  the  concentric  isothermal  picture  of  the  piire  heat  con- 
duct ionl^  so  that  almost  static  heat- conduct ion  relations  may,  here  be 
supposed.   An  estimate  of  the  Te  effect  made  under  this  assumption 
at  small  Gr  likewise  gives  good  agreement  with  the  measurements  (see 
following  section  b).   The  following  fundamental  difference  is  never- 
theless to  be  observed.   For  initially  given  cylinder  temperature  and 
heat  transfer,  the  solution  of  the  heat-conduction  equation  has  as  the 
potential  function  at  infinity  a  singular  point  (negatively  infinite 
temperature).   That  is,  even  at  a  very  large  distance  from  the  cylinder 
the  room  temperature  continues  to  decrease,  whereas  actually  in  the 
case  of  free  flow  the  room  temperature  at  large  distance  from  the  body 
very  soon  becomes  constant. 

(b)  Heat  transfer.  -  Dependence  on  Gr.   For  the  dependence  of  Nu 
on  Gr,  the  solution  of  the  differential  equations  under  the  familiar 
assumptions  of  the  Prandtl  boundary-layer  theory  gives 

Nu  -  Grl/4=  (6) 

as  is  also  approximately  shown  by  the  tests  of  Wamsler  and  Koch  (fig.  10), 
At  very  large  Gr  with  at  least  partial  turbulent  flow,  the  heat  trans- 
fer must  rise  more  strongly  on  account  of  the  increased  mixing.   With 
increasing  turbulent  mixing  the  dependence  of  the  heat-transfer  coef- 
ficient on  the  position  must  also  decrease.   The  assumption  that  a  is 
independent  of  the  position  (or  on  the  cylinder  diam.),  as  is  obtained 
from  the  tests  of  Griffiths  and  Davis  (ref.  17)  for  a  vertical  plate, 
gives  on  the  basis  of  dimensional  considerations 


1/3  (7) 


Nu  -  Gr 

Above  Gr  =  3.5^10  ,  the  proportionality  of  the  Nu  niunber  to  the 
fourth  root  of  Gr  corresponding  to  the  upstream  travel  of  the  point 
of  transition  to  turbulence  must  gradually  become  a  proportionality  to 
the  cube  root.   Experiments  of  W.  King  (ref.  18)  on  vertical  cylinders, 


1  7 
'See  photographs  by  R.  B.  Kennard;  Bur.  Stand.  J.  Res.,  vol.  8 

(1932),  p.  787.   Report  on  this  by  M.  Jakob:  Forschg.  Ing.-Wes.,  vol.  4 

(1933),  p.  45. 
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plates,  and  blocks  at  Gr  up  to  10   confirm  the  cube-root  law  as 
already  pointed  out  by  M.  Jakob  and  W.  Linke  (ref.  19).   For  very  small 
Gr,  for  which  the  amount  of  heat  transferred  by  convection  is  always 
below  that  transferred  by  conduction,  the  Wu  number  in  accordance 
with  its  definition  must  gradually  become  constant,  as  may  be  seen  from 
the  discussion  of  data  below  Gr  =  10"   (fig.  lO). 

Dependence  on  Pr.   The  dependence  of  Nu  on  Pr  for 
Gr  =  constant  and  Te  =  constant  may  be  understood  by  assuming  a 
change  in  X  with  all  the  remaining  magnitudes  kept  constant.   The 
increase  in  Pr,  through  a  decrease  in  X,  gives  in  the  region  of 
predominant  heat  conduction  (very  small  Gr)  a  proportional  decrease 
in  the  quantity  of  heat  transferred,  that  is,  Nu  =  constant.   In  the 
region  of  large  Gr,  a  decrease  in   X  (on  account  of  the  smaller  con- 
duction component  as  compared  with  the  convection  component)  gives  a 
relatively  smaller  decrease  in  a,  that  is,  an  increase  in  Nu.   In 
agreement,  investigations  by  Davis  (refs.  3  and  20 )  on  liquids  of  very 
different  Pr  (Pr  =  0.74  to  7940)  for  small  Gr  give  the  limiting 
value  of  Nu  which  is  independent  of  Gr  and  also  of  Pr.   For  medium 
and  large  values  of  Gr,  an  increase  of  Nu  with  Pr  was  obtained. 

Dependence  on  Te.   The  dependence  of  Nu  on  Te  for 
Gr  =  constant  and  Pr  =  constant  for  diatomic  gases  is  determined  by 
the  temperature  dependence  on  X,  v,  a,  and  also  by  the  temperature 
chosen  in  computing  the  characteristics  Gr  and  Nu  (Pr  is  independ- 
ent of  the  temperature)  for  the  fluid  properties   X  and  v  (P  =  l/Too 
is  constant  in  the  entire  field,  ref.  l);  the  wall  temperat\ire  t^  or 

T„  was  chosen  herein, 
w 

An  increase  of  Te  =  (T^/T^)-1  is  considered  to  be  due  to  a  de- 
crease in  T   with  T„  constant  so  that  the  values  of  X  and  v 

CO  W 

employed  for  the  computation  of  Nu  and  Gr  remain  unchanged.   The 
associated  increase  of  p  and  ®  is  compensated  by  a  decrease  in  the 
gravity  field  g   so  that   gP®,  that  is,  the  lift  acceleration  and  Gr 
remain  constant.   Hence  colder  outer  layers  with  smaller  X  and  v 
but  larger  p  now  take  part  in  the  heat  transfer.   In  the  region  of 
very  small  Gr  with  predominating  heat  conduction  X,  the  magnitude 
a,  and  hence  Nu,  must  therefore  decrease  with  increasing  Te.   This 
is  also  shown  by  the  data  discussed  previously  between  Gr  =  10"'^  and 
10  (see  fig.  10).   For  the  region  of  very  large  Gr  with  smaller  vis- 
cosity and  heat -conduction  effect  as  compared  with  the  turbulent  momen- 
tum and  heat  mixing,  it  also  follows-*-"  from  the  momentum  equation. 


1  R 

Under  the  reasonable  assumption  that  with  the  change  of  Te  the 
profiles  of  the  velocity  and  temperature  to  a  first  approximation  under- 
go affine  variations.   The  'standard  velocity  W  may  thus  be  taken  for 
example  as  the  maximum  value  of  the  tangential  velocity  at  the  equa- 
torial region  of  the  cylinder. 
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■because  of  the  constancy  of  the  lift  acceleration,  that  the  standard 
velocity  W  is  constant.   The  participation  of  colder  layers  with  in- 
creased value  of  p  in  the  heat  transfer  here  means  an  increase  in  the 
heat-transfer  coefficient  a,  which  is  proportional  to  cpW,  and  there- 
fore in  Nu. 

The  decrease  of  Nu  with  increasing  Te  for  very  small  Gr 
contrasts,  therefore,  with  an  increase  for  very  large  Gr.   The  effect 
of  Te  for  small  Gr  values  must  therefore,  in  agreement  with  the 
test  results  of  Langmuir,  first  decrease  with  increasing  Gr  values 
and  finally,  somewhere  in  the  range  of  medium  Gr,  must  practically 
vanish,  as  is  evidently  the  case  according  to  the  experiments  "between 
Gr  =  10^  and  lo'^  (fig.  10). 

For  the  range  of  very  small  Gr,  the  order  of  magnitude  of  the 
dependence  of  Nu  on  Te  under  the  approximating  assumption  of  purely 
static  heat  conduction  may  be  estimated  when  a  mean  conductivity  is 
used. 


(8) 


the  following  expression  is  obtained  for  the  Nu  numbers  corresponding 
to  two  different  Te  values  Te-^  and  Te2;  for  equal  wall  temperature 
T   =  T  • 


Nu(Te^)   X^(Tei) 

NuCTeg)  -^niCTeg)  ^^^ 


The  value  Te-]_  =0  is  chosen,  that  is,  vanishing  temperature  difference 
or  T^  =11;  and  Teg  =  1,  that  is,  T^  /T^  =  2.   For  the  temperature 

dependence  of   X,  the  previously  (ref.  1)  given  exponential  form  is 


dl9 


use 


with  n  =  0.738^° 


X/Xq  =  (t/Tq)"  (10) 


Xq  is  the  value  of  X  for  the  arbitrarily  chosen  temperatiore 


^0- 


Determined  from  the  previously  given  exponent  for  a  of  1.7  38 
neglecting  the  temperature  dependence  on  c. 
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From 


^m(Te2)  =  t! 


^2     2 


X  dT 


with  To  =  T   =  T^   in  equation  (lO)  and  by  integration: 


^0  -  "-M 


w 


\n(Te2)  =  ^TTT 


Mt^J   t^2 


wi> 


2    ^2     2 


T« 


"W- 


n  +  !■ 


Since  ^ni(Te-|_)  =  X(T^  ),  the  following  expression  results  after  sub- 
stitution of  the  numerical  values: 


X,(Te2) 


m^ 


'^^(Tiiy 


1.738 


-© 


1^1.738] 


=  0.81 


(11) 


The  decrease  of  Nu  for  Gr  =  constant  with  increase  from  Te  =  0 
to  Te  =  1  is  thus  theoretically  obtained  as  19  percent.   The  pre- 
viously discussed  tests  of  the  four  investigators  between  Gr  =  10~* 
and  10  gave  values  between  13  percent  and  26  percent  with  a  mean 
value  of  22  percent  (see  fig.  5).   This  good  agreement  justifies  the 
assumption  of  approximately  static  conducting  conditions  at  small  Gr. 


II.   THEORETICAL  SOLUTION  OF  THE  BOUNDARY  EQUATIONS  FOR  THE 

HORIZONTAL  CYLINDER  (CASE  OF  STEADY  MOTION) 

1.  Abstract 

The  large  number  of  experimental  investigations  of  heat  transfer 
in  free  convection  correspond  to  only  a  single  physically  satisfactory 
theoretical  solution  of  the  differential  equations,  namely  the  solution 
for  the  vertical  plate  (two-dimensional  steady-flow  case)  for  moderate 
temperature  differences,  which  was  given  by  Sclimidt  and  Beckmann 
(ref.  21)  with  the  help  of  Pohlhausen.   The  solution  is  based  on  the 
approximation  of  the  Prandtl  boundary -layer  theory  (laminar  flow  in  a 
layer  which  is  thin  as  compared  with  the  distance  from  the  lower  edge; 
velocity  normal  to  the  wall  small  as  compared  with  that  in  the  direction 
of  the  principal  flow  along  the  wall).   The  results  of  the  theory  as 
regards  the  velocity  and  temperature  fields  and  therefore  the  heat 
transfer  agree  very  well  with  the  corresponding  measurements. 
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For  the  case  of  the  horizontal  cylinder,  schlieren  photographs  in 
air  (see  ref.  22  and  figs.  23  and  25)  show  that,  for  Grashof  character- 
istics (Gr  =  d^gPS/vS)  of  about  10"^  and  above,  the  heating  of  the  air 
extends  to  only  a  relatively  thin  layer  around  the  cylinder,  a  fact 
which  permits  the  mathematical  simplification  of  the  problem.   On  the 
assumption  that  the  heat-transfer  and  the  flow  processes  are  restricted 
to  a  thin  film,  as  compared  with  the  cylinder  diameter,  with  laminar 
flow  (boundary-layer  assumption)  the  differential  equations  (ref.  23) 
for  moderate  temperature  differences  can  be  solved  approximately  for 
the  velocity  and  temperature  fields;  and  therefore  the  heat  transfer 
can  be  computed  which,  in  the  range  of  Gr  10  to  3x10  ,  within  which 
the  initial  assumptions  are  satisfied,  is  in  good  agreement  with 
experiment. 


2.  Setting  up  of  Differential  Equations 

To  start,  the  hydrodynamic  differential  equations  are  used  in  a 
form  previously  given  (ref.  l)   which  shows  a  "lift  term",  charac- 
terizing the  free  convection,  that  arises  from  the  combining  of  the 
gravity  term  and  the  hydrostatic  component  of  the  pressure  drop,  so 
that  only  the  gradient  of  the  dynamic  pressure  p*  remains  as  the 
pressure  force.   This  dynamic  pressure  is  equal  to  zero  where  there 
is  no  motion  and  no  temperature  difference  with  respect  to  the  medium 
at  a  large  distance.   As  was  indicated  in  detail  in  reference  1  (p. 
429),  the  density  and  other  characteristics  of  the  substance  in  the 
entire  field  may  be  considered  as  constant  (case  of  small  temperature 
difference)  and  the  equations  will  continue  to  describe  the  free- 
convection  problem.   The  equations  in  vector  form  are  then^S 

w°  grad  w  =  -  1/p  grad  p*  -  v  rot  rot  w  -  g|30 


div  w  =  0 
w°  grad  0  =  a  A 


(12) 


J 


with  the  boundary  conditions:  w  =  0,  0  =  ®  on  the  cylinder  surface;  and 
w  =  0,  0  =  0,  p*  =  0  at  infinity.   (w,  velocity;  0,  temperature  differ- 
ence with  respect  to  temperature  at  infinity;  p*,  dynamic  pressure; 
p,  density;  v,  kinematic  viscosity;  P,  expansion  coefficient;  a,  temper- 
ature conductivity;  g,  gravity  acceleration  vector;  g,  its  absolute 
value;  r,  cylinder  radius. ) 


In  equation  (3c)  in  reference  1  a  minus  sign  mistakenly  appears 
on  the  right  side. 

In  place  of  Aw  the  invariant  form  is  written  grad  div  w  - 

rot  rot  w. 
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For  treating  the  two-dimensional  flow  about  the  cylinder,  these 
equations  are  conveniently  written  in  arc  length  coordinates  with  the 
arc  length  s  on  the  cylinder  perimeter  (taken  in  clockwise  direction 
with  the  lower  stagnation  point  as  origin)  and  the  normal  distance  n 
from  the  cylinder  surface,  and  the  velocity  components  u  and  v 
parallel  and  normal  to  the  wall,  respectively.   There  is  then  obtained, 
without  neglecting  any  terms^'^: 


Su      Su     U'V        1  Sp 


» 


s 


u  -r—  +  V  ,-—  +  =  _  —  ^  +  gB0  sin  —  + 

r  +  n   OS     on   r  +  n     posr  +  n  r 

1-1    e • l/ e   1-6 

1  §^     1   Su      r^     S^u      u         2r    Sv 
[Sn^   r  +  n  Bn    (r  +  n)^  Ss2    (j-  +  n)^    (r  +  n)^  3s 
e2  i/e2    i/e  1        i        2  •  e 


r     Sv     Sv  u2  1  Bp*  „^     s 

u  r—  +  V  ■T—  -  =  -  —  3:—  -  gpe  cos  —  + 

r  +  nds     on  r  +  n  Pon           r 

1-6    e  .1  1 


L^    |i  +  |l+^L_.o  (13) 

+  n  OS   on   r  +  n  ^ 


r 

lie 


he         he        I  b^e  i      he  r2    h^e 


u  •5—  +  V  T^  =  a<  — -  +  -r—  +  ^   >  • 

r  +  n   5¥     hn         K  2   r  +  n  5^   (r  +  n)2  ^^2  / 

1-1    e-l/e  e2  i/e2        i/e  1 

The  approximations  made  according  to  the  boundary-layer  theory  are  now 
introduced.   The  cylinder  radius  r  and  hence  also  s   (with  the  ex- 
ception of  the  lower  stagnation  point)  and  u  are  assumed  to  be  of  the 
order  of  magnitude  1;  n  is  assumed  to  be  of  the  order  of  magnitude  of 


■^he  hydrodynamic  equations  agree,  except  for  the  dynamic  pres- 
sure and  lift  terms,  with  the  equations  (lO),  (ll),  and  (l2)  for 
r  =  constant  of  W.  Tollmien:  in  Wien-Harms,  Handb.  d.  Exp.  Phys.,  IV, 
1,  p.  248. 
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the  boundary-layer  thickness   e  «1^  so  that  n  +  r  may  be  re- 
placed by  r.   The  estimate  of  the  order  of  magnitude  of  the  individual 
terms  according  to  the  method  given  by  Prandtl  and  his  coworkers  then 
gives  the  results  written  underneath  the  individual  terms.   In  order 
that  the  friction  and  inertia  forces  may  be  of  equal  effect,  it  is 
necessary  that   v  (and  correspondingly  a)  be  of  the  order  of  magnitude 
of   e^.   In  order  that  the  lift  forces  for  the  motion  in  the  direction 
of  the  main  stream  be  of  significance  as  compared  with  the  inertia  and 
friction  forces,  that  is,  in  order  that  'free  convection'  exist  at  all, 
it  is  necessary  that  gPS  ■  1.   From  the  two  momentum  equations  it 
then  follows  that  the  pressure  terms  are  at  most  «=  1.   Since  p*  is 
zero  at  the  outer  edge  of  the  boundary  layer,  it  follows  from 

i  ^P*  «  1 
p  dn 

by  integration  from  the  outer  edge  of  the  boundary  layer  to  the  cylinder 
surface  over  the  length   e  that  ^ /q     in  the  boundary  layer  can  at 
most  be  of  the  order  of  magnitude   e,  that  is,  so  that  l/p  Sp*/5s  =  e. 


Since  in  what  follows  only  the  terms  of  the  order  1  are  retained, 
the  following  result,  important  for  simplifying  the  mathematical  treat- 
ment, is  obtained:  The  tangential  drop  of  the  dynamic  pressure  is  to 
be  neglected  as  compared  with  the  tangential  lift,  inertia,  and  friction 
forces.   The  term  p*  alone  is  contained  in  the  equation  of  the  normal 
momentum.   Hence,  u,  v,  and  0  are  to  be  computed  from  the  equation  of 
the  tangential  momentum,  the  continuity  equation,  and  the  heat-transfer 
equation  (eq.  (l4)).   Thereafter  p*  may,  if  desired,  be  computed  from 
equation  (15)  of  the  normal  momentum.   Hence,  the  boundary-layer  equa- 
tions for  the  horizontal  cylinder  for  free  convection  are 


Su 
ds 


au 


=  V 


.2 

O  U       nn  ■ 

—  +  gpe  sm 


Su   Sv  _ 
Ss"   5n  ~ 


he         'be 

u  V-  +  V  V- 
os     on 


=  a 


an2 


A 


(14) 


; 


U'" 

r 


1  Sp«    „„     s 

V—  -  gPS  cos  — 

p  dn    '^'^      r 


(15) 
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A  corresponding  estimate  of  equations  (l3),  first  rendered  non- 
dimensional,  gives  the  two  corresponding  conditions: 

1     2   J  Gr    ^   .    _,     1  .    . 

_.,2^nd  —  .  1,  I.e.  Gr-  -  (b) 

For  a  maximum   g,  in  the  sense  of  the  approximations,  of  0.1  (approxi- 
mation to  about  10  percent),  this  gives  Re  =  100  and  Gr  =  10^  as 
lower  limits  of  the  computations.   The  approximation  is  therefore  better 
the  larger  Gr,  that  is,  the  smaller  the  boundary-layer  thickness  e . 
The  upper  limit  is  then  attained  if  the  boundary  layer  becomes  turbu- 
lent, which  (according  to  part  IV)  is  the  case  for  Gr  =  3.5X10®  at 
the  upper  stagnation  point,  for  Gr  =  SxlO^  at  the  equator.   The  range 
of  applicability  of  the  laminar -boundary-layer  computations  therefore 
lies  between  Gr  =  10'^  and  Gr  =  10^. 


3.  Reduction  to  Ordinary  Differential  Equations 
The  usual  introduction  of  the  stream  function  \|f 

eliminates  the  continuity  equation.   By  introducing  nondimensional 
magnitudes  with  Gr'  =  r^gps/v^  "by  means  of  the  transformations^; 


s 
X  =  — 

r     "    r 


y=IlGr'l/4      Y=iGr'-l/4     ^=^      (l?) 


there  remain  of  the  five  constants  (v,  gP,  a;  0,  r)  of  differential 
equations  (14)  only  two  in  the  nondimensional  combination  of  the  Prandtl 
characteristic  Pr  =  v/a. 

5y  dx  dy  "  5x  ^y2  ~  ^^3 

(18) 

Sy  Sx   Sx  Sy   Pr  Sy2 

Equation  (l7)  is  not  an  arbitrary  transformation  which  carries 
equation  (14)  over  into  equation  (l8),  but  a  definitely  determined 
transformation  with  the  aid  of  which  it  is  possible  to  eliminate  Gr 


Gr'   referred  to  r  is  more  convenient  for  the  computation  than 
Gr  referred  to  d.   In  the  final  result,  the  more  usual  Gr  is  again 
introduced. 
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so  that  it  no  longer  occurs  in  equation  (l8).   The  solution  of  equa- 
tion (18 )  therefore  likewise  no  longer  contains  Gr.   It  is  therefore 
unnecessary  to  know  the  solution  of  equation  (iS)  if  it  is  of  interest 
to  know  only  how  the  solution  (e.g.,  velocity,  temperature,  boundary- 
layer  thickness,  heat-transfer  coefficient)  of  the  initial  equation  (14) 
depends  on  the  Gr  numher.   This  dependence  of  Gr  is  already  com- 
pletely represented  by  the  transformation  (l7),  whereas  equation  (iS) 
contains  the  further  dependence  on  Pr  and  on  the  space  coordinates. 
In  other  words,  the  solution  of  equation  (iS),  obtained  for  a  definite 
Pr   value24a  gives,  by  means  of  equation  (l7),  at  the  same  time,  the 
velocity  and  temperature  fields  for  all  Gr  values.   On  the  other  hand, 
the  transformation  (17)  without  the  solution  of  equation  (iS)  already 
gives  important  general  information  on  the  flow  condition  and  the  heat 
transfer  for  free  convection  with  respect  to  the  dependence  on  Gr. 
From  equation  (17)  there  follows,  if  the  following  expression  is  set  up, 


that 


Il  =  yGr'-l/4   ^  =  uGr'V2   ^^vGr'VS   £|9=|Ig^,1/4 
r   ■'  V  V  9  on   oy 

(20) 

That  is,  the  relative  boundary -layer  thicknesses  decrease  with  Gr"  /  ; 
the  nondimensional  velocities  increase  with  Gr  /  ;  and  the  non- 
dimensional  temperature  drop  increases  with  GrV'^,  from  which  there 
follows  directly  the  well-known  l/4-power  law  of  heat  transfer  in  free 
convection: 


V^  (21) 


Nu  ■  Gr 

These  theoretical  results  agree  with  those  for  the  rectangular 
plate  and  are  not  restricted  to  the  latter  and  the  horizontal  cylinder 
alone,  but  depend  on  more  general  assumptions  (ref.  19).   They  are 
valid,  in  general,  for  all  two-dimensional  cases  of  free  flow  with  a 
boundary-layer  character  (ref.  23),  that  is,  where  the  thickness  of  the 
heated  laminar  layer  is  small  compared  with  the  distance  from  the  inci- 
dence edge  and  the  radius  of  curvature  of  the  wall,  and  the  latter  does 
not  change  dis continuously  in  the  flow  direction.   The  existence  of 
sharp  edges  about  which  the  flow  occurs  is  therefore  excluded.   In  this 
more  general  case,  there  are  on  the  right  side  of  equation  (l3)  addi- 
tional terms  which  are  determined  by  the  radius  of  curvature  and  which 


It  is  here  assumed  as  self  evident  from  physical  considerations 
that  one  and  one  only  such  solution  exists.   On  the  mathematical  side 
of  these  existence  and  uniqueness  proofs,  see  the  discussion  in  refer- 
ence 1  (p.  428). 
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vary  with  the  arc  length  (Sr/Ss)  but  which  drop  out  in  the  estimate 
made  under  the  initial  assumptions.   Furthermore,  in  place  of  the  sine 
and  cosine  of  the  lift  term  there  appear  functions  which  determine  the 
direction  of  the  surface  element  with  respect  to  the  horizontal  and 
vertical  but  which  similarly  do  not  change  anything  in  the  transfor- 
mations (17)  and  the  conclusions  from  equations  (20)  and  (2l)  based  on 
these  transformations. 

In  agreement  with  these  theoretical  conclusions  equation  (21)  has 
already  been  confirmed  for  vertical  plates  (ref.  21),  horizontal  cyl- 
inders^S,  and  cylindrical  layers  (ref.  24),  and  even  for  the  heat 
transfer  from  the  two  sides  of  a  square  horizontal  thick  plate  with 
edges  (ref.  25)  and  for  the  heat  transfer  with  evaporation  at  a  verti- 
cal cylinder  in  water  and  carbon  tetrachloride  (ref.  26)  for  which  the 
previously  mentioned  assumptions  are  not  all  satisfied. 

A  reduction  of  the  partial  differential  system  (I8)  to  systems  of 
ordinary  differential  equations  is  attempted  by  setting 


for  example, 


q  =  y  g(x)  (22a) 

;x,y)  -  p(q)f(x)  (22b) 

T(x,y)  =  t(q)  (22c) 


U  =  ^  =  f(x)g(x)p'(q)  (22d) 


based  on  the  notion  that  all  the  profiles  of  the  stream  function  ^(y), 
tangential  velocity  U(y),  and  temperature  T(y)  for  the  different 
azimuths  x  of  the  cylinder  are  produced  by  affine  distortions  by 
means  of  the  azimuth  functions  f(x)  and  g(x)  from  a  single  "base 
profile"  of  the  stream  function  p(q)  or  the  velocity  p'(q)  and  one 
base  profile  of  the  temperature  t(q)  (q  =  normal  coordinate  of  the 
base  profile).   Substitution  of  equation  (22)  in  equation  (I8)  gives 
(primes  of  p  and  t  denote  differentiation  with  respect  to  q;  primes 
of  f  and  g  denote  differentiation  with  respect  to  x): 

P'2(f2gg'  +  ff'g2)  -  pp"ff'g2  =  p"'fg3  +  t  sin  x 

(23) 
l/Pr  gt"  +  f'pt'  =  0 


'^^See  Part  I,  section  5,  figure  10. 
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A  separation  into  functions  of  the  two  independent  variables  x 
and  q,  that  is,  a  successful  application  of  the  assumed  expression 
(eq.  (22))  for  solving  the  equations,  is  possible  only  if  the  two  func- 
tions f (x)  and  g(x)  can  be  uniquely  determined  from  the  four  ordinary 
differential  equations 

f'(x)  =  ag(x)  (24a) 

f(x)f' (x)g2(x)  =  c  sin  X  (24b) 

f2(x)g(x)g' (x)  =b  sin  X  (24c) 

f(x)g3(x)  =  d  sin  X  (24d) 

(where  a,b,c,d  are  initially  undetermined  constants)  in  spite  of  the 
redundancy  of  the  equations^^.   The  possibility  of  such  unique  determi- 
nation, which  is  the  necessary  assumption  for  the  applicability  of 
equation  (22),  will  be  discussed  in  more  detail  after  equation  (29)  is 
considered.   If  such  determination  is  possible,  there  remain  for  p(q) 
and  t(q)  two  ordinary  differential  equations: 

(b  +  c)p<2(q)  -  cp(q)p"(q)  =  dp'"  (q)  +  t(q)     1 

(25) 
t"(q)  +  a  Pr  p(q)f  (q)  =  0  J 

In  section  4,  equation  (24)  will  be  solved  and  in  section  5  equation  (25) 
will  be  discussed. 

4.  Determination  of  the  Azimuth  Functions  F(x)  and  G(x) 

(a)  Setting  up  of  the  two  differential  equations  for  F(x).  -  The 
boundary  values  of  f(x)  and  g(x)  which  are  required  for  solving  equa- 
tion (24)  are  determined  from  physical  considerations.  From 

St   d.t(q)   ,  V 
57  =  ^  e(x) 

it  follows  that 

g(x  =  0)  =  go  ^  0  (26a) 


The  author  is  indebted  to  Dr.  A.  Waumann,  Leipzig,  for  several 
mathematical  suggestions. 
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Otherwise,  the  temperature  of  the  cylinder  surface  would  also  he  the 
temperature  in  the  field  along  the  entire  normal  at  the  lower  stagna- 
tion point,  which  cannot  "be  the  case.   From  equation  (24a)  there  is 
then  ohtained 

f'(x  =  0)  =  agQ  (26h) 

Since  the  tangential  velocity  at  the  lower  stagnation  point  must,  for 
reasons  of  symmetry,  "be  different,  there  follows  from  equation  (22d): 

f(x  =  0)  =  0  (26c) 

The  system  of  equation  (24)  with  the  "boundary  conditions  (26a),  (26"b), 
and  (26c)  thus  contains  five  free  availa"ble  constants  a,b,c,d,  and  gQ, 

which  are  determined  in  subsequent  computation. 
By  introducing  the  normed  functions^ ' : 

F(x)  =  ^M      G(x)  =  SM  (27) 

^§0  So 

equation  (24)  may  be  transformed  to: 

F'  =  G  (28a) 

f2  G  G'  =  J^   sin  x  (28b) 

F  F'  G^  =  ^   sin  x  (28c) 

a  e 

F  g2  =  — ^  sin  X  (28d) 

ago 

with  the  boundary  conditions 

F(x  =  0)  =  0      F'(x  =  0)  =  1      G(x  =  0)  =  1       (29) 

A  necessary  assumption  for  the  applicability  of  the  expressions  of 
equation  (22)  for  solving  the  equations  is  that  the  two  functions  F 
and  G  can  be  uniquely  determined  from  the  four  equations  (28a)  to 
(28d) .   This  would  be  the  case  if  it  could  be  proven  that  the  four 


27 

By  normed  is  meant  fully  determined  also  with  respect  to  a 


factor. 
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equations  are  not  independent  of  one  another,  "but  that  two  equations 
follow  from  the  other  two.   This  proof  is  partly  possible  only  inso- 
far as  one  equation  (for  suitable  choice  of  the  constants)  is  a  con- 
sequence of  the  remaining  equations.   Thus,  elimination  of  G  with 
the  aid  of  equation  (28a)  gives 

f2  pi  p"  ^  J^        sij^  X  (30a) 

F  F'3  =  J^   sin  X  (30b) 

^     d 

F  F'-J  =  J  sin  X  (30c) 

ago* 

Equation  (30c)  is  thus  a  consequence  of  equation  (30b)  in  the  case  that 
the  following  first  condition  for  the  constants  is  satisfied: 


d  =  c/a  (31) 


Setting  for  briefness 


c/o2 


/e?g^^   =  a  (32a) 

b/a^g^^:  =  p  (32t) 

there  then  remain  for  the  determination  of  F  the  two  equations 

F  F'3  =  a  sin  X  "       (33l) 

p2  p,  p'.  =  p  sin  X  (33II) 

with  freely  available  values  of  a  and  p.   As  a  necessary  assumption 
for  the  applicability  of  equation  (22),  there  remains  the  further  con- 
dition that  the  two  equations  I  and  II  possess  two  equal  (or  at  least 
approximately  equal)  solutions. 

This  proof  will  be  obtained  by  solving  equations  I  and  II  by  a 
power-series  development  above  and  below  the  equator  x  =  jt/2.  For 
this  purpose,  with  cp  =  x  -  n/2,  I  and  II  are  transformed  into 

F((p)F'2(<p)  =  a  cos  <p  (34l) 

f2((p)F' (<p)F"((p)  =  p  cos  (p  (34II) 

where  the  primes  now  denote  differentiation  with  respect  to  (p. 
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(b)  Solution  of  differential  equation  (l)  for  F.  -  The  value  of 
a,  by  the  boundary  values  of  equation  (29),  is  first  determined  as 

a  =  1  (32c) 

as  is  seen  from  the  first  term  of  the  power-series  development  of 
equation  (331 )  about  the  point  x  =  0.   From  the  expression  for  solving 
equation  (34l): 

F^(»)  =  aQ(l  +  C;]_  ^  +  cg  (p2  +  .  .  .)  (35) 

there  follows 

F^(»  =  0)  =  aQ      F-'-'(<p=  0)  =  aQ  c^^  (36a) 

and  by  substituting  in  equation  (341) 

aQ4  c-L^  =  a  =  1  (36b) 

Substitution  of  expression  (35)  in  differential  equation  (34l) 
gives  the  unknown  coefficients,  which  are  taken  up  to  the  fifth,  in- 
clusive, and  which  can  be  successively  determined  in  terms  of  c-.  : 

^2  =  -  I  =1^     =3  =  -  ^  ^1  +  ^  =1^    =4  =  ^  ^1^  -  ^  ^1^ 

"(37) 

1            5       ,      13      p;  '^    ' 

Cr   =  —  Ct   -  Ct  -^  +  Cn"3 

5     360   1   324   1    216   1 

The  boundary  condition  F-'-(9  =  -:n:/2)  =  0  according  to  equation  (29) 
gives  for  the  value  of  c-^  the  solution  of  an  algebraic  equation  of 

the  fifth  degree  which  is  solved  by  trial 

c-^  =   0.581  (38) 

Of  the  five  possible  real  roots  it  is  necessary  to  choose  that  one  for 
which  no  further  zero  of  F^   is  obtained  between  (p  =  0  and   <p  =  -rt/2. 
The  coefficients  then  assume,  according  to  equation  (37),  the  values 

c-L  =  +  0.581     ^2  =  -  0.05626     Cj  =  -  0.01412   "j 

(39) 
C4  =  -  0.00165     ^5  =  -  0.00066  l 

According  to  equations  (36a),  and  (36b),  there  is  obtained 

F^{<f  =  0)    =  a-Q   =  1.504  F^'((^=  0)    =  aQ  c-^   =  0.874  (40) 
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and  the  numerical  solution  of  F-'-((p)  is  thereby  completely  givBn  (see 
table  III). 

It  may  be  remarked    that  for  equation  (331)  a  solution  was  also 
obtained  by  a  power-series  expansion  at  the  point  x  =  0  in  powers  of 
sin  (x/2)  which,  because  of  the  fact  that  the  curvature  of  F  at  that 
point  is  zero,  converges  with  particular  rapidity: 

F*(x)  =  2  sin  I  (l  .  3^  sin2  |  .  ^  sin4  |  ^  .  .  .)   (4l) 

This  solution  assumes  at  the  equator  the  values 

F*(x  =  3t/2)  =  1.500     F*'(x  =  rt/2)  =  0.874  (e) 

which  are  in  good  agreement  with  the  values  of  equation  (40 )  of  the 
series  expansion  at  the  equator.   In  this  way,  the  sufficient  conver- 
gence of  fI(<P)  with  the  five  computed  terms  is  assured,  and  a  good 
check  is  obtained  against  computation  errors.   The  solution  F*  will 
not  be  used  in  what  follows;  this  is  because  of  the  singular  behavior 
of  equation  (33Il)  at  the  point  x  =  0,  and  the  identity  proof  must 
therefore  be  carried  out  by  the  series  expansion  at  the  equator. 

(c)   Solution  of  differential  equation  II  for  F.  -  For  equa- 
tion (34II),  corresponding  to  equation  (35),  the  following  expression 
is  assumed: 

F^^(,p)    =  bQ(l  +  d-^   +   dg*^  +  .  .  .)  (42) 

Substitution  in  the  differential  equation  gives  the  following  values 
for  the  constants:  bg  and  d.-^     are  arbitrary 

d  -    P 


2  ~  2bo4d^ 


2 


(43) 


'^l  S  =  -  I  ^^1  ^^2  +  ^2) 

where,  on  account  of  the  length  of  the  expressions,  d^  and  d^  are 
not  given  as  functions  of  d-[_  and  dg.   The  function  F    thus  depends 


?7a  •• 

Acknowledgment  is  made  to  Prof.  W.  Muller,  Aachen,  for  the 
suggestion  that  the  solution  of  equation  (331)  may  be  reduced  to  two 
quadratures  by  separation  of  the  variables.   One  of  the  quadratures  is 
not  in  closed  form,  however,  but  can  be  carried  out  only  by  a  series 
expansion  so  that  this  method  would  not  be  essentially  different  from 
the  one  employed. 
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on  the  as  yet  freely  available  values  of  Tdq,  d-^,    and  |3.   The  solu- 
tions F   and  F    are  equal  if  the  corresponding  coefficients  of  the 
series  are  equal.   Therefore,  the  as  yet  undetermined  coefficients  are 


equated  as  follows ; 


The   condition 

d2   = 

=  =2 

Tdq  =  aQ 

gives 
3 

21^0*^1 

di  =  c^  (44) 


L 
or  with  equations  (36"b)  and  (44) 

|3  =  -  l/3  a  =  -  1/3  (45) 

All  the  coefficients  may  now  be  expressed  in  terms  of  d-|_: 

^2  =  -  I  <il'    ^3  =  +  fe  ^l'    ^4  =  +  72  ^l'  -  T2   ^1^ 

Comparison  of  the  coefficients  of  the  two  series  for  F   and  F 
gives 

bo  -  ao  =  0    d^  -  c^  =  0    dg  -  cg  =  0   d3  -  C3  =  ^  c^ 

^4  -  ^4  =  -  2^  =^2    d5  -  C5  =  -^  3^^  ci  +  ^  0^3       (47) 

This  result  states  that  the  solutions  of  the  two  different  differen- 
tial equations  (34l)  and  (34Il)  agree  at  the  equator  in  the  value  of  the 
function,  in  the  tangent,  and  in  the  curvature,  but  not  in  the  higher 
derivatives  in  which  they  differ  fundamentally  to  a  small  degree.   F-'- 
and  F^-'-  are  approximately,  but  not  rigorously,  equal.   The  system  (24) 
is,  strictly  speaking,  overdetermined  and  has  only  approximately  unique 
solution.   This  solution  of  expression  (22)  is  not  rigorously,  but 
only  approximately,  valid.   Because  of  the  agreement  of  F-^  and  F 
at  the  equator  up  to  and  including  the  curvature,  however,  these  differ- 
ences become  more  appreciable  in  the  neighborhood  of  the  lower  and 
upper  stagnation  points  where  the  initially  assumed  boundary-layer 
assumptions  are  not  satisfied.   The  method  that  has  here  been  employed 
thus  gives  an  approximate  solution  of  the  boundary-layer  equations  with- 
in the  accuracy  limits  detennined  by  the  physical  approximations  of  the 
boundary-layer  assumptions. 
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The  numerical  computation  of  equations  (44)  and  (46)  with  Bq 
and  C3_  according  to  equations  (40)  and  (39)  gives: 

ho  =  1.504    d-j_  =  0.581    "^2   "   ~   0-05626    ^^  =  +   0.01816 

d^^  =  -0.00321   d.5  =  +  0.00144  (48) 

Table  III  shows  the  values  of  the  functions  fI((^)  and  Fll((p)  computed 
according  to  equations  (39),  (40),  and  (48),  with  sufficient  agreement 
between  x  =  30  ,  150  and  even  at  165   (graphically  interpolated  for 
F   ).   The  negative  value  of  F    for  x  =  0  has  no  significance 
since,  on  account  of  the  singular  behavior  of  equation  (34II)  for 
F  =  0,  the  series  (42)  beyond  this  singular  point  no  longer  represents 
the  solution  of  the  differential  equation.   For  the  further  equations, 
F   is  employed  because  only  this  function  possesses  the  required 
property  of  vanishing  at  the  lower  stagnation  point. 

(d)  Determination  of  G.  -  The  term  G  is  determined  according 
to  equation  (28d)  where,  according  to  equation  (49),   d  =  agQ'i. 

Table  III  shows  the  computed  values  based  on  F-"-. 


5.  Determination  of  Basic  Profiles  p(q)  and  t(q) 

The  constants  b,  c,  d  which  occur  in  equation  (25)  can  now,  with 
the  aid  of  equations  (31),  (32a),  to  (32c),  and  (45)  be  expressed  in 
terms  of  a  and  gQi 

b  =  -  1/3  a2  gQ4    c  =  a2  gQ4    d  =  a  gQ*         (49) 

Equation  (25)  then  becomes 

2/3  a2  gQ^  p'2  _  a2  gQ^  p  p"  =  a  gQ*  p"'  +  t   | 

)  (50) 

t"  +  a  Pr  p  t'  =  0  J 

with  the  boundary  conditions: 

q=0:p=0     p'=0     t=l 

q  =  •:  p'  =  0     t  =  0  (g) 

Equation  (50)  for  a  =  1  and  gg  =  1  may  now  be  numerically  integrated. 

The  following  values  were  chosen  for  the  as  yet  free  remaining  constants 
a  and  gQ 
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a  =  3      gQ  =  3-1/4  (5^) 


equation  (50)  becomes 

p"'  +3pp"-2p'2  +  t=0 
t"  +  3  Pr  p  t'  =0 


(52) 


This  system  with  the  same  boundary  conditions  and  for  Pr  =  0.733 
(diatomic  gases)  has  already  previously  been  numerically  integrated  by 
E.  Pohlhausen  for  the  vertical  plate  so  that,  because  of  the  introduc- 
tion of  the  normed  functions  F  and  G  and  the  consequent  availability 
of  the  constants,  a  new  numerical  integration  is  unnecessary.   The 
stream  function  p(q)  and  temperature  function  t(q)  correspond  to  the 
functions   C(C)  and  0(C)  of  Pohlhausen.   Their  values  are  given  in 
reference  21  in  tabular  form  (ref.  4,  p.  187).   Graphical  representa- 
tions of  the  velocity  profile  p'(q)  and  of  the  temperature  profile 
t(q)  are  given  in  figures  13  to  16  (continuous  curves). 


6.  Complete  Solution  for  Velocity  and  Temperature  Fields 

With  the  values  of  equations  (51)  for  a  and  gQ,  the  azimuth 

functions  f(x)  and  g(x)  according  to  equation  (27)  are  now  uniquely 
determined.   For 

a  gQ  =  3^4  =  2.280      gQ  =  3-^^^  =  0.760         (53) 

there  is  obtained 

f(x)  =  2.280  F(x)      g(x)  =  0.760  G(x)  (54) 

Table  TV  shows  the  values  of  f(x),  g(x),  and  f(x)-g(x)  computed  from 
F   and  G  of  table  III.   Figure  17  shows  these  values  graphically. 

With  the  values  f (x)  and  g(x)  from  table  IV,  the  values  p(q) 
and  t(q)  from  the  table  of  Pohlhausen  and  the  transformations  (l6), 
(17),  and  (22)  there  is  obtained  the  complete  solution  of  the  boundary - 
layer  equations  (14)  for  the  velocity  (u  and  v)  and  temperature 
fields  (0)  obtained  without  any  empirical  value.   Collecting  results 
yields  (see  also  eq.  (l")  and  (20)  for  simultaneous  passage  from  Gr' 
to  Gr  (Gr  =  8Gr') : 
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u(s,n)  =  - 


v(s,n) 
0(s,n> 
with  X 


V  Grl/2 


> 


V  GrV^  1^  .  ,  f    .    f(x)g'(x)     ,,  , 
-7-^|f'(^)P(^)  ^   g(x)     '^P'(q) 

et(q)      Ih^"^^)  =|^s(x)t'(q) 


s   ^     n  GrV^  ,    . 

-   and  q  = t-t—  g(x) 

r     ^   r  qI/'^ 


(55a- 
)   55d) 


y 


For  the  determination  of  the  normal  velocity  v  it  is  required  to 
know  f ' (x)  and  g'(x).   The  term  f'(x)  is,  according  to  equation  (24), 
determined  by  means  of  g(x);  g' (x)  was  determined  by  graphical  differ- 
entiation of  g(x)  °.   Here  the  qualitative  description  of  the  normal 
velocity  obtained  from  the  streamline  picture  is  sufficient.   It  can  be 
shown  that  the  boundary  conditions  of  u,  v,  0  for  q  =  0  are  cor- 
rectly assumed.   If  the  assumption  n  <c r  which  was  made  with  respect 
to  the  solution  of  equation  (14)  is  applied  for  very  large  q,  there 
are  obtained  for  u  and  0  the  correct  value  zero;  but  for  v,  in 
contradiction  to  the  boundary  conditions  (12),  a  finite  value  is  ob- 
tained which  is  related  to  the  fact  that  the  cylinder  curvature  was 
neglected. 

The  boundary-layer  thickness,  tangential  velocity,  temperature, 
and  heat-transfer  coefficients  according  to  equation  (55)  behave  as 
follows.   With  regard  to  the  dependence  of  these  magnitudes  on  Gr, 
there  hold  first  of  all  the  considerations  of  section  3  in  connection 
with  equations  (20)  and  (2l).   The  dependence  on  the  cylinder  azimuth 
is  such  that  all  profiles  of  the  tangential  velocity  and  temperature 
are  obtained  by  affine  distortions  from  the  basic  profiles  p'(q)  and 
t(q). 

The  value  l/g(x)  (fig.  17)  represents  the  'extension'  of  the  tem- 
perature and  velocity  profiles  normal  to  the  surface  with  increasing 
azimuth,  that  is,  the  development  of  any  characteristic  distance  from 
the  wall,  for  example,  of  the  place  of  maximum  velocity  (q  =  0.95)  or 
of  the  boundary-layer  thickness  to  be  defined  later  (q  =  2.18).   In 
accordance  with  this,  the  boundary-layer  thickness  at  the  lower  stag- 
nation point  possesses  a  finite  approximately  constant  value  over  a 
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The  values  are  not  given  here. 
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large  azimuth  range  but  increases  at  first  gradually  then  more  and  more 
rapidly  and  attains  at  the  upper  stagnation  point  a  theoretically  in- 
finitely large  value  (upward  current  of  warm  air).  The  function  g(x), 
since  it  is  proportional  to  the  temperature  drop,  likewise  represents 
the  variation  of  the  local  heat-transfer  coefficient  with  the  cylinder 
azimuth.   At  the  lower  stagnation  point,  g(x)  is  constant  and  with 
constant  slope  connects  with  the  values  of  the  other  side  of  the  cyl- 
inder.  At  the  upper  stagnation  point  (x  =  it)  this  is  no  longer  the 
case;  g(x)  has  a  sharp  peak  which  is  associated  with  the  fact  that  the 
boundary-layer  assumptions  at  this  point  strongly  deviate  from  the 
actual  conditions  (normal  velocity  large  as  compared  with  the  tangen- 
tial velocity).   Under  actual  conditions  this  peak  is  balanced  out  by 
a  very  sharp  minimum  value  as  is  shown  by  the  value  of  Nu  along  the 
cylinder  perimeter  computed  from  the  schlieren  pictures  of  E.  Schmidt 
(ref.  22)   .   No  quantitative  comparison  of  this  experimental  curve 
of  g(x)  with  the  theoretical  curve  given  herein  can  be  given  because 
the  conditions  of  the  two-dimensional  problem  are  so  little  satisfied 
for  the  schlieren  pictures  (pipe  length  =  2  times  cylinder  diam. )  that 
the  computed  heat -transfer  coefficient  is  40  percent  greater  than  that 
of  the  tests  of  Koch  (ref.  22). 

The  development  of  the  tangential  velocity  with  the  azimuth,  that 
is,  the  increase  of  any  characterizing  velocity,  for  example  the  max- 
imum velocity,  is  described  by  f(x)'  g(x)  (see  fig.  17).   The  velocity 
increases  from  the  value  zero  at  the  lower  stagnation  point  (up  to  60° 
approximately  linearly)  to  a  maximum  at  about  128°  and,  on  approaching 
the  upper  stagnation  point  (convergence  of  the  flow  on  either  side) 
rapidly  decreases  to  zero.   This  physically  correct  behavior  of  the 
solution  at  the  upper  stagnation  point  is  not  introduced  as  a  boundary 
condition  but  is  obtained  as  a  necessary  consequence  of  the  theory. 

An  over-all  picture  of  the  velocity  and  temperature  fields  is 
given  by  the  streamlines  and  isotherms  (figs.  11  and  12).   They  are 
computed  from  the  equations 


tllz^lll  =   f(x)Gr'l/4  p  [£  Gr'lA  g(x)]  =  const 

®        Lr        ^^  ^J 

for  various  values  of  the  constants30.   The  scale  of  the  representation 


(55e,f) 
const 


^^See  figure  13  for  maximum  Gr  =  16x10^. 

The  streamline  picture  for  azimuth  intervals  of  15°  and  in  addi- 
tion for  5°,  10°,  170°,  175°;  the  isotherm  picture  for  azimuth  intervals 
of  30°  and  in  addition  for  165°,  170°,  175°. 
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holds  for  Gr  =  10  ;  for  Gr  =  10   the  distances  from  the  cylinder 
surface  should  be  reduced  by  a  factor  of  3.16  and  for  Gr  =  10^,  by  10. 
At  the  lower  and  equatorial  neighborhood  of  the  cylinder  up  to  an  azi- 
muth of  about  105°,  the  flow  near  the  wall  is  directed  toward  the  sur- 
face, while  above  it  is  directed  away  from  the  surface  upwards.   In 
this  respect  the  theoretical  solution  for  the  cylinder  differs  funda- 
mentally from  that  for  the  vertical  plate,  which  gives  only  a  flow 
toward  the  plate.   The  isotherms  at  the  lower  part  form  almost  con- 
centric circles.   At  the  upper  stagnation  point  (upward  warm  current) 
the  isotherms  theoretically  go  toward  infinity,  whereas  actually  (corre- 
sponding to  the  gradual  dissolution  and  spreading  out  of  the  warm  up- 
streaming  air)  they  close  at  a  relatively  large  distance  over  the 
cylinder.   The  isotherm  0.53  gives  at  the  same  time  the  place  of  the 
maximum  tangential  velocity. 


7.  Heat-Transfer  Law 

From  the  expression  for  the  quantity  of  heat  dQ  passing  from 
an  element  of  area  df, 


dQ 


=  X(-^  [n  =  0])  df  (56) 


the  value  of  the  temperature  gradient  at  the  wall  obtained  from 
equation  (55c) 

^  (n  =  0)  =|Gr>lA  g(x)t>o  (57) 

(t'Q  =  t'(q  =  O)  =  -0.508),  and  the  defining  equation  for  locally 
variable  heat-transfer  coefficient  a(x) 

a(.)=i§  (58) 

there   is  obtained 

2^^=  g(x)(-to')  Gr'lA  (59) 

From  this  there  is  obtained  for  the  dependence  of  the  nondimensional 
local  heat-transfer  coefficient  on  the  cylinder  azimuth  x  and  the 
Gr  number  (Gr^  =  8Gr')^°^: 


Since  several  formulas  of  this  section  will  be  used  later  in 
Part  III,  and  in  order  to  avoid  misunderstandings,  indices  (e.g.,  d,h,H) 
for  the  characteristics  Gr  and  Nu  are  written  which  indicate  the 
lengths  with  which  these  magnitudes  are  formed. 


32  NACA  TM  1366 

2:1^^  =   NU(i(x)  =0.604  g(x)  Gr^lA  (60) 

see  figure  18. 

To  compute  the  mean  heat-transfer  coefficient  averaged  over  the 
cylinder  perimeter, 

a(x)  dx  (61) 


there  is  required  the  mean  value 

(62) 

which,  "by  planimetering  g(x)  (fig.  17),  is  found  to  "be 

g  =  0.616  (63) 

It  may  also  be  computed  hy  the  following  relation  from  equation  (24a) 
with  a  =  3  according  to  equation  (51) 

3  rt  g  =  f (jt)  (64) 

which  gives  g  ~   0.620.   However,  the  value  0.616  is  preferable  because 
this  value  takes  into  account  the  total  curve  g(x)  and  does  not  depend 
on  the  end  point  of  the  series  expansion  f(x).   From  equation  (SO)  the 
following  is  then  obtained  for  the  mean  heat-transfer  coefficient  as  a 
function  of  Gr  for  free  flow  at  a  horizontal  cylinder 


2i^  =  Nu^  =  0.372  Gr^V^  (65) 


A  comparison  of  this  theoretical  heat-transfer  coefficient  for  diatomic 
gases  (Pr  =  0.74)  with  test  results  is  given  in  Part  I,  section  5;  see 
also  figure  10. 


8,  Comparison  of  Theoretical  Velocity  and  Temperature 

Fields  with  Available  Measurements 

The  comparison  of  the  computed  velocity  and  temperature  fields  is 
carried  out  with  the  measurements  of  Jodlbauer  (ref.  16)  which  are  the 


2r  =  9   cm, 

"t^    —       yv.c,    f 

t   =  18.1^, 

2r  =  5   cm. 

t^  =  104.6°, 

t   =  18.1°, 
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31 
only  ones  thus  far  available"^  .   Of  the  fields  measured  for  six  differ- 
ent  Gr   values,  the  following  ones  were  computed: 

Gr  =  3.76xlo6  (figs.  13,  14) 
Gr  =  6.54x105  (figs.  15,  16) 

The  best  agreement  with  theory  is  to  be  expected  for  the  maximum  Gr. 
On  the  basis  of  the  theory  it  would  be  desirable  to  have  measurements 
for  the  largest  possible  laminar  Gr,  that  is  (see  Part  IV)  values  of 
about  3.5x10  ,  which  could  be  realized  in  air  at  20°  C  at  a  surface 
temperature  of  100°  C  with  cylinders  of  42  centimeters  diameter. 

The  comparison  with  the  theory  was  carried  out  in  such  manner  that 
all  velocity  profiles  u(s,n)  and  temperature  profiles  0(s,n)  corre- 
sponding to  a  value  of  Gr  were,  with  the  aid  of  equation  (55),  re- 
computed for  the  theoretical  basic  profiles  p'(q)  and  t(q)  in  the 
representation  of  which  they  must  all  coincide. 

In  the  case  of  the  velocity  profiles  (figs.  13  and  14),  there  is 
a  regular  deviation  in  that  the  measured  velocities  are  greater  than 
the  computed  ones,  the  deviation  from  the  theory  in  the  maximum  velocity 
amounting  on  the  average  to  22  percent  for  the  smaller  and  17  percent 
for  the  larger  value  of  Gr.   A  part  of  the  deviations  occurring  for 
large  distances  from  the  wall  beyond  the  maximum  value  should  be 
ascribed  to  the  uncertainties  in  the  difficult  measurements  of  such 
small  air  velocities,  as  can  be  seen  from  the  following  discrepancies. 
The  measured  velocities  at  the  azimuth  60°  for  Gr  =  3.76x10^   lie 
considerably  above  the  theoretical  curve,  whereas  for  Gr  =  6.54x10^ 
they  lie  on  the  theoretical  curve.   On  the  other  hand,  the  velocities 
measured  at  x  =  30°  for  Gr  =  6.54x10^  show  strong  deviations  upward, 
whereas  for  Gr  =  3.76x10°  they  are  in  agreement  with  the -theory.   The 
measured  distance  from  the  wall  of  the  maximum  value  of  the  velocity  is 
everywhere  in  very  good  agreement  with  the  theory. 

The  measured  temperature  profiles  show  very  good  agreement  with 
the  theory  even  for  the  azimuth  of  165°.   A  small  regular  deviation  is 
noted  for  all  azimuths  and  Gr  numbers  for  medium  distances  from  the 
wall.   The  small  deviations  at  small  distance  from  the  wall  can  be 
ascribed  partly  to  the  fact  that  the  condition  for  the  isothermal  sur- 
face was  not  quite  satisfied  in  the  tests.   The  measured  temperature 
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drops  at  the  wall  for  Gr  =  6.54x10^  are  on  the  average  3  percent 
smaller  than  the  values  given  by  the  theory,  while  for  Gr  =  3.76  xLO^ 
they  are  1  percent  higher-^^. 

It  is  likewise  of  interest  to  compare  the  solution  based  on  the 
boundary-layer  theory  for  the  vertical  plate  with  the  corresponding 
test  results.   Reference  21  (figs.  20  to  23)  also  shows  the  velocity 
fields  to  have  stronger  deviations  than  the  temperature  fields.   For 
the  small  plate  with  smaller  deviations,  the  measured  velocities  also 
lie  somewhat  higher  than  the  theoretical  values;  but  for  the  larger 
plate  with  larger  deviations  the  contrary  is  true.   For  the  temperature 
fields,  particularly  for  the  large  plate,  no  regular  deviation  can  be 
established  between  theory  and  experiment.   The  deviations  occurring 
for  the  plate  are  clearly  smaller  in  comparison  with  those  for  the 
cylinder.   This  is  to  be  ascribed  to  the  fact  that  in  the  case  of  the 
plate  only  the  boundary-layer  terms  of  the  order  e2  <<  i  (rectilinear 
flow)  were  neglected,  whereas  for  the  case  of  the  cylinder  there  were, 
in  addition,  neglected  terms  of  the  order  el  <<  1  that  were  related 
to  the  curvature . 


III.   HYDRODYNAMIC  AND  THERMAL  COMPARISON  BETWEEN  VERTICAL 

PLATE  AND  HORIZONTAL  CYLINDER  FOR  FREE  CONVECTION 

AND  FOR  PLATE  IN  PARALLEL  FLOW 

1.  Abstract 

Now  that  the  theoretical  boundary -layer  solutions  for  the  vertical 
plate  and  the  horizontal  cylinder  for  free  convection  have  been  obtained, 
it  is  of  advantage  to  compare  these  two  standard  bodies  of  two- 
dimensional  free  flow  with  regard  to  the  shape  of  the  stream  and  the 
heat  transfer.   Similarly,  a  comparison  of  the  boundary-layer  develop- 
ment between  the  vertical  plate  for  free  convection  and  for  the  plate 
in  parallel  flow  (Blasius  solution  of  the  boundary-layer  equations)  is 
of  interest.   The  following  features  are  characteristic  for  the  three 
types  of  flow. 


The  fact  that  the  values  of  the  heat-transfer  coefficient  ob- 
tained from  the  field  measurements  of  Jodlbauer  through  integration  of 
the  temperature  gradients  at  the  wall  over  the  cylinder  perimeter  lie, 
on  the  average,  4  percent  higher  than  the  theoretical  values  (see 
fig.  lO)  is  due  to  the  deviation  between  the  theoretical  and  the  actual 
conditions  in  the  neighborhood  of  the  upper  stagnation  point. 
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The  boundary -layer  thickness  of  the  plate  in  parallel  flow  in- 
creases with  the  square  root  of  the  distance  from  the  incidence  edge; 
that  of  the  vertical  plate  in  free  convection  increases,  however,  as 
the  fourth  root  of  the  height  above  the  lower  edge  (cf.  eqs.  (69)  and 
(71)).   Both  start  with  the  thickness  zero,  which  for  the  plate  in 
parallel  flow  results  in  infinitely  large  velocity  gradients,  and  for 
the  free  plate  results  in  an  infinitely  large  temperature  drop  and  an 
infinitely  large  local  heat-transfer  coefficient  at  the  lower  edge. 
In  the  case  of  the  cylinder,  on  the  contrary,  the  boundary  layer  starts 
at  the  lower  stagnation  point  of  the  horizontal  cylinder  with  a  finite 
thickness  and  thereafter  with  finite  velocity  and  temperature  gradients 
and  a  finite  local  heat-transfer  coefficient  (fig.  17).   It  increases 
according  to  a  complex  law  (l/g(x),  see  eqs.  (70)  and  (82))  and  at  the 
upper  stagnation  point  reaches  a  theoretically  infinite  thickness,  with 
vanishing  velocity  and  temperature  drops  normal  to  the  surface  (upward 
stream  of  warm  air)  and  with  a  vanishing  heat -transfer  coefficient. 

In  spite  of  these  fundamental  differences  it  is  possible  to  set 
up  a  number  of  relations  between  the  three  cases.   These  show  on  the 
one  hand  a  close  hydrodynamic  and  thermal  kinship  which  is  due  to  the 
boundary- layer  character  of  the  differential  equations  underlying  the 
theoretical  solutions.   On  the  other  hand,  they  give  practical  view- 
points with  regard  to  the  application  of  rectangular  plates  or  hori- 
zontal cylinders  for  the  heat  transfer.   For  this  purpose  it  is  first 
of  all  necessary  to  represent  the  Reynolds  numbers  of  the  characteristic 
length,  of  the  boundary-layer  thickness,  and  of  the  nondimensional  local 
and  average  heat-transfer  coefficients  as  functions  of  the  Grashof 
numbers  of  the  plate  and  cylinder,  which  in  turn  requires  the  introduc- 
tion of  a  boundary-layer  thickness  (flow  discharge  thickness)  for  the 
velocity  profile  of  the  free  convection.   A  knowledge  of  the  dependence 
of  Re  of  the  boundary -layer  thickness  on  Gr  of  the  plate,  as  well 
as  of  the  azimuth  and  Gr  of  the  cylinder,  is  in  addition  required 
for  the  later  evaluation  of  the  tests  as  regards  the  occurrence  of 
turbulence  (Part  IV)  for  free  convection  at  the  plate  and  the  cylinder. 


2.  Boiondary-Layer  Thickness  for  Free  Convection 

It  is  first  necessary  to  determine  the  magnitude  to  be  associated 
with  the  boundary -layer  thickness.   For  the  velocity  profile  u(y)  of 
the  plate  in  parallel  flow  (U  =  maximum  velocity),  Prandtl  and  his 
coworkers,  as  is  known,  introduced  the  displacement  thickness  5*, 
defined  as 

6*.  U=      (U  =  u)  dy  (66) 

Jo 

which  is  a  measure  of  the  decrease  in  the  flow  discharge  as  a  result 
of  the  friction. 
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For  the  velocity  profile  of  the  free  convection  u(y),  a  flow 
discharge  thickness  5  is  introduced  which  is  hydrodynamically  equiv- 
alent to  the  displacement  thickness,  being  like  the  latter  character- 
istic of  the  development  of  the  friction  layer,  which  in  this  case 
takes  up  the  entire  flow: 


U  = 


u  dy  (67) 


or  expressed  in  words:   Through  the  flow  discharge  thickness  6  would 
flow  the  same  fluid  mass  with  the  maximum  velocity  as  actually  flows 
with  the  total  stream. 

In  the  nondimensional  velocity  profile  p'(q)  (q  =  nondimensional 
distance  from  the  wall,  p  =  stream  function)  which  occurs  in  the  theo- 
retical solution,  the  discharge  flow  thickness  is  denoted  "by  qg,  the 
maximum  velocity  hy  P™a.x*   From  equation  (67)  there  is  then  obtained: 

q^P'  ^  =      P'(q)dq  =  P(-)  -  P(0) 


6  max 


10 


and  with  the  numerical  values: 

p'    =  0.275 
max 

p(c»)  =  0.60^^ 

p(0)    =  0 
the   value   of      qg      is   obtained  as 

qg   =  2.18  (68) 

For  orientation  purposes^^^  it  is  assumed  that  the  velocity  maximum 
lies  at  q  =  0.95  and  the  point  of  inflection  of  the  profile  at 
q  =  1.85. 

With  the  values  of  q^,  for  the  plate  and  cylinder  from  equa- 
tions (72)  and  (82),  respectively,  there  is  obtained  for  the  vertical 
plate,  where  h  is  the  distance  from  the  lower  edge: 


•Extrapolation  from  last  computed  value  p(6.0)  =  0.5928. 
See  continuous  curves  in  figures  13  and  15. 


NACA  TM  1366  37 

I^Gr^l/4  ^2.18  (69) 

and  for  the  horizontal  cylinder  of  diameter  d  at  azimuth  x 

I  g(x)2V4  Gr//^  =  2.18  (70) 

3,  Nondimensional  Representation  of  Theory  for  Vertical  Plate 

(a)  Reynolds  numbers.  -  From  the  theoretical  solution  for  the 
vertical  plate  given  hy  Schmidt  and  Pohlhausen  (ref.  21)  the  following 
expression  is  obtained  according  to  their  equations  (23),  (24),  and 
(31)'^^,  since  c^  =  Grj^/4h^,  for  the  velocity  u  and  the  distance 

from  the  wall  y: 

n=^V'  (q)  Gr^l/S      ^  =  q  V^  Gr^'lA  (7l) 

For  the  flow  discharge  thickness  y  =  5  and  the  maximum  velocity 
u  =  U,  there  is  then  obtained 

U=XPiaxGrhV2  |  =  q^  ^2   Gr^-l/^  (72) 

or 

^  -  2  V^  Piax  %  Gr^^/^  ^  =   2  P '  ^x  ^r.V^  (73) 

and  using  the  abbreviations  R^  =  U-h/'^   and  Rg  =  U-&/v   and  substi- 
tuting the  numerical  values  for  Pj^j^x  ^^^     ^ 


Rg  =  1.695  Grjj 


1/4  (74) 

R^  =  0.550  Grj^V^  (75) 

Elimination  of  Grj^  from  the  last  two  equations  gives  the  following 

expression  as  the  relation  between  Re  of  the  flow  thickness  and  the 
height  of  the  plate  for  a  vertical  plate  in  free  flow: 


The  distance  from  the  lower  edge  is  denoted  by  h  instead  of  x, 
the  nondimensional  distance  from  the  wall  ^   by  q,  the  nondimensional 
stream  function  t,     by  p. 


38  NACA  TM  1366 


Rg  =  2.29  'R^f'^  (76) 

(b)  Nusselt  characteristics.  -  For  the  local  heat-transfer  coef- 
ficient a(h)  at  the  height  h  ^aliove  the  lower  edge,  the  following 
equation  holds  (ref.  2l)*: 

1/4 


1  fGrvAV* 
a(h)  =X(-t,.)i(^-/j  (77) 


(t_'  =  -0.508  is  the  value  of  the  derivative  of  the  nondimensional 

temperature  function  t(q)  for  q  =  O)  or,  in  nondimensional  form,  the 
numerical  values  may  be  substituted: 

5Mlh  ^  0.359   Gr^lA  (78) 

The  variation  of  the  heat-transfer  coefficient  along  the  variable 
height  h  of  a  plate  of  total  height  H  is  most  conveniently  obtained 
by  introducing  a  local  heat-transfer  coefficient  which  has  been  made 
nondimensional  by  dividing  through  by  H 

^   =  NuH(h)  .  0.359  GrnV^  (|)"'/'  (79) 

The  heat-transfer  coefficient  (fig.  18)  at  the  lower  edge  of  the  plate 
with  zero  boundary-layer  thickness  is  theoretically  infinitely  large 
and  then  continuously  drops  with  the  reciprocal  of  the  fourth  root  of 
the  height. 


The  mean  heat-transfer  coefficient  corresponding  to  a  plate  of 
height  H 

1 


am  =  H 


H 

a(h)dh  (j) 

0 


is  obtained  from  equation  (77)  as 

^  =  |-  #M-to')|GrHlA  (80) 

or,  nondimensionalized  and  with  n\imerical  values  substituted  (fig.  18), 


as 


^  =  Nujj  =  0.479   GrgV^  (81 ) 


*  Extrapolation. 


MCA  TM  1366  .  39 


4.  Nondlmensional  Representation  of  Theory  for  Horizontal  Cylinder 

(a)  Reynolds  numbers.  -  The  theoretical  solution  given  previously 
for  the  tangential  velocity  u  (eq.  (55a))  and  the  distance  n  from 
the  wall  (eq.  (55d))  for  the  azimuth  x  of  a  horizontal  cylinder, 
with  the  values  n  =  6  and  the  maximum  velocity  u  =  U,  may  he  trans- 
fonned  to 

U  =  I   f  (x)g(x)pi,,  Gr//2      I  =  -^  Gr,-l/4        (82) 

For  the  upper  stagnation  point  x  =  it,  "because  g(jt)  =  0  and 
f(jr)  =  5.84,  the  boundary -layer  thickness  becomes  infinite  and  the 
tangential  velocity  zero.   The  quantity  flowing  through  U6  and  there- 
fore the  Re  of  the  boundary  layer  U&/v,  however,  remain  finite.   To 
introduce  Re  for  the  characteristic  length  results  in  no  simplication 
for  the  cylinder  because,  in  addition  to  the  characteristic  length,  the 
cylinder  radius  occurs  as  an  additional  characteristic  length  of  the 
system.   From  equation  (82), 


U5/v  =  f(x)p'  ^  q„  Gr^V^  (83) 


-max 


Substituting  the  numerical  values  and  using  Gr,^  result  in  the  follow- 
ing expression  for  Re  of  the  boundary  layer  at  azimuth  x: 

7(q{x)    =   0.357  f(x)  Gr^V^  (84) 

and  for  Re  of  the  boundary  layer  at  the  upper  stagnation  point 
(with  f(rt)  =  5.84): 

R5  =  2.08  Gr^jV^  (85) 

where  Rg  is  written  for  Rg(n).   Although  the  theoretical  solution 

for  X  =  jt  is  no  longer  valid,  equation  (85)  because  of  the  small 
change  of  f (x)  gives  the  correct  relations  in  the  region  of  the  upper 
stagnation  point  (up  to  about  165°),  so  that  for  the  sake  of  simplic- 
ity the  formula  will  be  used  in  the  following  discussion. 

(b)  Nusselt  characteristics.  -  Figure  18  shows  the  variation  of 
the  nondlmensional  local  heat-transfer  coefficient  according  to  equa- 
tion (6O)  over  the  developed  semicircumference  of  the  cylinder.   In 
the  region  of  the  lower  stagnation  point,  it  has  a  finite,  practically 
constant  value;  and  at  the  upper  stagnation  point,  it  drops  rapidly  to 
zero  in  accordance  with  the  upstreaming  warm  air  at  that  point.   In  the 
same  figure  the  mean  heat-transfer  coefficient  Nu^  for  a  cylinder 

is  plotted  according  to  equation  (65). 
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5.  Hydrodynamic  Comparison 

(a)  Betveen  vertical  plate  and  horizontal  cylinder  in  free  con- 
vection. -  Although  the  boundary- layer  growth  for  a  vertical  plate  and 
a  horizontal  cylinder,  as  already  remarked  (Part  III,  sec.  l)  and  as 
seen  from  equations  (72)  and  (82),  differ  fundamentally  in  character, 
it  is  nevertheless  of  interest  to  compare  the  flow  condition  (always 
characterized  in  what  follows  by  Rg)  at  the  upper  edge  of  the  plate 

and  in  the  region  of  the  upper  stagnation  point  of  the  cylinder.   The 
following  three  questions  will  here  be  considered: 

1.  The  behavior  of  the  flow  at  the  upper  edge  of  the  plate  and  at 
the  upper  stagnation  point  of  the  cylinder  for  the  case  that  the  height 
of  the  plate  is  equal  to  the  cylinder  diameter. 

2.  The  behavior  of  the  flow  at  the  upper  edge  of  the  plate  and  at 
the  upper  stagnation  point  of  the  cylinder  for  the  case  that  the  height 
of  the  plate  is  equal  to  the  developed  semicircumference  of  the  cyl- 
inder, that  is,  the  characteristic  lengths  are  equal. 

3.  The  relation  of  the  height  of  a  plate  whose  flow  condition  at 
the  upper  edge  is  equal  to  that  of  the  upper  stagnation  point  of  a 
cylinder  to  the  diameter  of  the  cylinder. 

Before  each  of  these  questions  is  considered,  the  ratio  of  the 
boundary-layer  Re  numbers  at  the  upper  stagnation  point  of  the  cyl- 
inder R5(z)  and  at  the  upper  edge  of  the  plate  R5(P)  are  written, 

with  the  aid  of  equations  (74)  and  (85): 

R5(Z)/R5(P)  =  1.226  Qv^/^/Gt^/^  (86) 

where  R5(z)  and  Gr^^  refer  to  the  cylinder,  and  R5(P)  and  Qtv-^ 

refer  to  the  plate.   The  coefficient  1.226  has  the  following  signifi- 
cance (cf.  eqs.  (85)  and  (83)  for  the  cylinder,  eqs.  (74)  and  (73)  for 
the  plate ) : 


2.08 


2~^/^Pmax-%-f(^)     fU) 


'•'''-^^"   23/2p.    .   -^  ^^'^ 

■^max  ^5 

1.  For  equal  Gr  of  plate  and  cylinder,  Gr^^  =  Grj^,  equation  (86) 
gives 

R5(Z)  =  1.226.R5(P)  (88) 
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that  is,  for  equal  Gr  of  plate  and  cylinder,  for  example,  for  a  plate 
of  the  height  of  the  cylinder  diameter  h  =  d  under  otherwise  equal 
conditions  (with  regard  to  temperatures  and  materials),  Re  at  the 
upper  stagnation  point  of  the  cylinder  is  22  percent  greater  than  that 
of  the  upper  edge  of  the  plate  (fig.  19).   This  result  states  nothing 
about  the  hydrodynamic  relation  between  the  two,  since  the  comparison 
refers  to  different  characteristic  lengths.   It  is  useful  only  for  the 
rapid  comparison  of  plate  and  cylinder  for  the  usual  values  of  Gr. 
It  is  then  known  that  for  equal  Gr  the  cylinder  possesses  the  more 
developed  flow,  so  that  the  flow  at  the  cylinder  may  be  turbulent  while 
that  of  the  plate  is  still  laminar. 

2.  For  equal  characteristic  length  h  =  rt/2d  there  is  obtained 
from  equation  (86) 

R5(Z)/R5(p)  =  1.226  (2/rt)V^  =  0.874  (89) 

The  boundary-layer  Re  at  the  upper  stagnation  point  of  a  cylinder  is 
thus  13  percent  smaller  than  that  of  a  plate  of  the  height  of  the 
developed  semicircumference  of  the  cylinder  (fig.  19).   This  retar- 
dation of  the  boundary-layer  development  as  compared  with  the  vertical 
plate  is  due  to  the  fact  that  the  boundary  layer  at  the  lower  stagnation 
point  of  the  cylinder  already  has  a  finite  thickness,  whereas  in  the 
case  of  the  plate  the  thickness  must  increase  from  zero. 

3.  For  equal  boundary-layer  Re   values  for  the  plate  and  cylinder, 
equation  (86)  gives  » 

Gr^^Gr^  =  1.226"'^  =  0.441  (90) 

At  the  upper  stagnation  point  of  the  cylinder,  there  is  therefore  the 
same  flow  condition  as  at  the  upper  edge  of  the  plate  if  Gr  of  the 
cylinder  is  0.44  times  that  of  the  plate.   This  permits  a  rapid  con- 
version of  the  flow  data  for  the  plate  into  those  for  the  cylinder  and 
vice  versa.   If,  for  example,  it  is  known  from  heat-transfer  experiments 
on  a  plate  that  the  departure  from  the  laminar  GrV^  law  occurs  at 
(Grj^)j^^  =  8^108,  then  it  can  be  concluded  without  experiment  from  equa- 
tion (90)  that  this  must  be  the  case  for  the  horizontal  cylinder  for 
('^•'"d^kr  ~  3.5x10  .   If  the  plate  and  cylinder  are  under  otherwise  equal 
conditions,  it  is  possible  from  equation  (90)  to  derive  a  simple  rela- 
tion for  the  height  of  a  plate  hg^  which  is  hydrodynamically  equivalent 

to  a  cylinder: 

d3:hg^3  =  0.441      h-^:d  =  1.311  (91) 
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In  the  neighborhood  of  the  upper  stagnation  point  of  a  horizontal  cyl- 
inder, the  same  flov  condition  prevails  (Rg)  as  under  other\^^ise  equal 

conditions  at  a  vertical  plate  of  height  1.31  times  that  of  the  cylin- 
der diameter  (fig.  19).   As  was  to  be  concluded  from  the  answers  to 
questions  1  and  2,  this  height  must  lie  between  the  cylinder  diameter 
and  the  developed  semicircumference. 

The  answer  to  the  three  questions  on  the  flow  condition  at  a 
cylinder  and  at  plates  of  various  heights  is  sho\rni  in  figure  19  (Rg 

of  the  cylinder  set  equal  to  l). 

(b)  Hydrodynamic  comparison  between  a  vertical  plate  in  free  con- 
vection  and  a  plate  in  a  parallel  flow.  -  The  relation  (76)  between  Re 
of  the  boundary-layer  thickness  and  that  of  the  height  of  the  plate 
(characteristic  length)  R^^  for  the  vertical  plate  in  free  convection 

is  strongly  analogous  to  the  corresponding  relation  (92)  between  Re 
of  the  boundary- layer  thickness  Rg*(referred  to  the  displacement 

thickness  5*)  and  the  characteristic  length  R^^  for  a  plate  in  a 

stream  parallel  to  its  plane  (Blasius  solution  of  the  boundary -layer 
equations) : 


R? 


2.29  •  Rh"^/^  (76)  II  Rg*  =  1.73  ^^/'^  (92) 


In  spite  of  the  already  mentioned  entirely  different  rate  of  growth  of 
the  boundary -layer  thickness,  in  the  case  of  the  parallel  flow  as  the 
square  root  of  the  distance  from  the  incident  edge  and  in  the  case  of 
the  free  convection  as  the  fourth  root  of  distance  from  the  lower  edge, 
the  development  of  the  Reynolds  number  of  the  boundary-layer  thickness 
as  a  function  of  the  Reynolds  number  of  the  characteristic  length  in 
the  case  of  these  two  very  different  types  of  flow  follows  the  same 
exponential  law  and  with  approximately  equal  magnitude. 


6.  Thermal  Comparison  between  Vertical  Plate  and  Horizontal  Cylinder 

(a)  Mean  heat-transfer  coefficients.  -  For  equal  Gr  of  plate  and 
cylinder  (Grg  =  Gr^),  there  is  obtained  from  equations  (8l)  and  (65)  the 

following  ratio  of  the  nondimensional  heat-transfer  coefficients  of  the 
horizontal  cylinder  Nu^^  and  of  the  vertical  plate  Nug 

NUj^/Nug  =  0.777  (93) 

Under  otherwise  equal  conditions,  therefore,  the  mean  heat-transfer 
coefficient  of  a  cylinder  is  only  78  percent  of  that  of  a  plate  of  the 
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height  of  the  cylinder  diameter  (fig.  18).  In  agreement  vith  this,  the 
comparison  given  by  M.  Jakob  and  W.  Linke  (ref .  19)  of  convection  tests 
on  vertical  plates  and  horizontal  cylinders  in  a  single  Wu-Gr  diagram 
(with  H  and  d  denoting  the  length  dimension)  gives  a  lower  position 
of  the  cylinder  tests.  The  modified  determination  of  the  critical  num- 
ber, which  follows  from  these  results,  will  be  discussed  later  (Part  IV, 
sec.  5). 

For  equal  mean  heat-transfer  coefficient  a^  of  plate  and  cylinder 

under  otherwise  equal  conditions,  there  is  obtained  from  equations  (8I) 
and  (65) 

0.479  H"V^  =  0.372  d'V"^      H/d  =  2.76  (94) 

The  mean  heat-transfer  coefficient  of  a  cylinder  is  therefore,  under 
otherwise  equal  conditions,  equal  to  that  of  a  plate  of  height  2.76 
times  that  of  the  cylinder  diameter.   On  account  of  the  small  decrease 
of  the  local  heat-transfer  coefficient  at  large  plate  heights,  this 
factor  of  2.76  is  not  very  sharply  determined.   Thus  Jodlbauer  (ref.  29) 
found  it,  possible,  without  considering  the  previously  given  theoretical 
solution  for  the  cylinder,  to  represent  the  mean  heat-transfer  coef- 
ficients for  cylinders  from  tests  by  Koch  by  the  formulas  for  the  ver- 
tical plate  if  H  is  replaced  by  2d  (instead  of  theoretically  by 
2.76d).   As  a  matter  of  fact,  the  a^  values  thus  determined  for  the 

cylinder  lie  only  8.3  percent  above  the  values  given  by  the  theory 
according  to  equation  (65),  as  can  be  seen  by  replacing  H  by  2d  in 
equation  (8I) . 

(b)  Total  heat  transfer.  -  The  total  heat  transfer  Q(z)  of  a 
cylinder  of  diameter  d  along  the  entire  circiimference  and  the  total 
heat  transfer  Q(P)  along  the  two  sides  of  a  plate  of  the  height  H, 
under  otherwise  equal  conditions,  are  in  the  ratio 

Q(Z)  ^  am(z)rtd  ^  ^   Nua(z) 
QXpT  "  ani(P)2H  "  2  Nuh(P) 

and  from  equations  (8I)  and  (65), 

Q(Z)/Q(P)    =  1.221  Gr^l/VGrjjl/4  (95) 

where  the   significance  of  the  factor  1.221  is  as   follows    (cf.    eqs.    (Sl), 
(80),    (65),    and   (59)): 

rt  0.572        jt   2^/^g(-to')  5Tci  ,      ^ 

1-22^  =  2  07479   =  2  4^  =  IW^  ^''^ 

3     2      ^"^0  -' 
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A  comparison  of  this  factor  1.221  with  the  factor  1.226  occurring  in 
the  ratio  of  the  boundary-layer  Re  of  cylinder  and  plate  (eq.  (86)) 
according  to  equation  (87)  shows  that  both  are  identical,  since  accord- 
ing to  the  theoretical  solution,  the  following  relation  holds  as  an 
immediate  consequence  of  the  validity  of  the  heat-transfer  equation: 

3rtg  =  f(jt)  (97) 

The  small  difference  of  l/2  percent  is  explained  by  the  fact  that  (as 
mentioned  in  connection  with  eq.  (64))  g  was  computed  not  according 
to  equation  (97)  but  was  obtained  by  planimetering. 

For  equal  Gr^^  and  Grij  values  in  equations  (86)  and  (95),  the 
following  is  thus  obtained: 

Q(Z)/Q(P)  =  R5(Z)/R5(P)  (98) 

that  is,  the  total  heat  transfer  of  a  horizontal  cylinder  and  that  of 
the  two  sides  of  a  vertical  plate  under  other^^rise  equal  conditions  are 
in  the  same  ratio  as  the  boundary- layer  Re   values  at  the  upper  stag- 
nation point  of  the  cylinder  and  at  the  upper  edge  of  the  plate, 
respectively. 

On  the  basis  of  relation  (98)  it  is  now  possible,  without  further 
computation,  to  answer  the  three  questions  on  the  total  heat  transfer 
of  a  cylinder  and  plates  of  different  height,  analogous  to  the  questions 
referring  to  the  hydrodynamic  comparison  (see  fig.  19;  the  heat  transfer 
of  the  cylinder  has  been  set  equal  to  l). 

Corresponding  to  equation  (88), 

Q(Z)  =  1.22  Q(P)  (99) 

The  total  heat  transfer  of  a  cylinder  is  thus  22  percent  greater  than 
that  for  the  two  sides  of  a  plate  of  a  height  equal  to  the  cylinder 
diameter  under  otherwise  equal  conditions  (fig.  18). 

Corresponding  to  equation  (89), 

Q(Z)  =  0.87  Q(P)  (100) 

The  total  heat  transfer  of  a  cylinder  is  13  percent  smaller  than  that 
of  the  two  sides  of  a  plate  of  height  equal  to  the  developed  semi- 
circumference  of  the  cylinder  (equal  characteristic  length)  under  other- 
wise equal  conditions. 
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Corresponding  to  equation  (91)  for 

Q(Z)  =  Q(P) 

there  is  obtained  the  height  of  a  plate  H^  thermally  equivalent  to 
a  cylinder: 

1%  :  d  =  1.31  (101) 

The  total  heat  transfer  of  a  cylinder  under  otherwise  equal  conditions 
is  therefore  as  great  as  for  the  two  sides  of  a  vertical  plate  of 
height  1.31  times  that  of  the  cylinder  diameter.   In  summary  it  may 
therefore  be  stated  that  a  horizontal  cylinder  and  a  vertical  plate, 
which  transfers  heat  on  both  its  sides,  of  a  height  equal  1.31  times 
the  cylinder  diameter  are  equivalent  both  thermally  (with  reference  to 
the  total  heat  transfer)  and  hydrodynamically  (with  reference  to  the 
boundary-layer  development  at  the  upper  stagnation  point  and  at  the 
upper  edge,  respectively). 

IV.   THE  OCCUREIENCE  OF  TURBULENCE  IN  FREE  FLOW  ABOUT  A  HORIZONTAL 

CYLINDER  AND  ALONG  A  VERTICAL  PLATE 

1.  Introductory  Observations 

For  the  determination  of  the  transition  from  laminar  to  turbulent 
flow  in  the  case  of  free  convection,  tests  were  carried  out  on  a  verti- 
cal plate  and  a  horizontal  cylinder  of  sufficient  size.   Two  main 
reasons  underlay  the  investigation.   In  the  first  place,  it  was  desired 
to  determine  the  upper  limit  of  the  validity  of  the  laminar-  boundary- 
layer  theories  for  plate  and  cylinder  and  therefore  the  upper  limit  of 
applicability  of  the  heat-transfer  formulas  developed  from  these  theo- 
ries for  practical  application.   The  second  reason  was  of  a  more 
theoretical  nature.   The  numerous  turbulence  investigations,  both  ex-' 
perimental  and  theoretical,  undertaken  in  recent  decades  have  been 
concerned  almost  exclusively  either  with  the  flow  between  two  parallel 
walls  with  linear  velocity  distribution  (Couette  flow)  or  with  the  flow 
in  pipes,  channels,  about  rotating  cylinders,  along  plates,  cylinders, 
or  other  resistance  bodies  which  have  in  common  a  velocity  profile  which 
rises  uniformly  from  the  value  zero  at  the  wall  up  to  a  maximum  value. 
Differing  essentially  from  these  profiles  are  evidently  velocity  pro- 
files for  which  the  velocity  rises  from  the  value  zero  at  the  wall  to 
a  maximum  value  and  then,  after  passing  through  a  point  of  inflection, 
again  drops  to  the  value  zero  at  a  large  distance  from  the  wall.   Such 
profiles  are  typical  in  the  case  of  heat  transfer  by  free  convection. 
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Some  further  observations  on  the  setting  up  of  turbulence  in  free  con- 
vection will  be  made  later  in  connection  with  the  discussion  of  results 


(sec.  5), 


2.  Test  Procedure 


For  the  investigation,  the  schlieren  method  developed  by  E.  Schmidt 
was  applied,  which  as  an  optical  method  provides  the  possibility  of  an 
instantaneous  view  of  the  entire  field,  and  is  therefore  of  advantage 
for  determining  the  critical  number. 

The  schlieren  method  has  been  described  in  detail  by  E.  Schmidt 
(ref.  22)  so  that  only  the  essential  points  will  be  given  here.   It  is 
based  on  the  deflection  suffered  by  a  light  ray  in  passing  through  a 
field  with  density  stratification  (i.e.,  variation  of  the  index  of 
refraction)  normal  to  its  direction  of  propagation  toward  the  colder 
layer.   In  the  test,  a  parallel  light  beam  is  allowed  to  pass  tangen- 
tially  along  the  surface  of  the  heated  body.   The  rays  in  the  neighbor- 
hood of  the  wall,  because  of  the  maximum  temperature  drop  at  that  place, 
are  most  strongly  deflected  from  their  initial  direction;  those  farther 
away  from  the  wall  are  less  deflected,  while  those  rays  which  are  out- 
side the  temperature  field  undergo  no  deflection.   On  a  screen  which 
is  set  up  at  a  sufficient  distance  behind  the  heated  body,  the  rays 
near  the  wall  appear  the  farthest  from  the  wall.   The  following  picture 
is  obtained:  the  heated  body  throws  a  completely  dark  shadow  which,  as 
compared  with  that  of  the  cold  body,  is  increased  by  the  thickness  of 
the  temperature  boundary  layer.   This  shadow  is  bounded  by  a  first 
"interior  caustic  curve"  adjacent  to  which  is  formed  a  medium-bright 
region  (rays  which  have  passed  through  the  temperature  field),  which  in 
turn  is  bounded  toward  the  outside  by  a  second  "external  caustic  curve" 
originating  from  the  tangential  rays  in  the  neighborhood  of  the  heated 
surface.   Still  farther  toward  the  outside  is  the  uniformly  bright  field 
of  the  undisturbed  illumination. 

For  an  accurate  evaluation,  the  method  is  suitable  only  for  bodies 
having  one  of  its  length  dimensions  (two-dimensional  problem),  so  that 
the  effect  of  the  temperature  fields  at  the  ends  on  the  path  of  the 
light  ray  is  small  compared  with  the  deflection  undergone  in  passing 
through  the  distance  along  the  length.   On  the  other  hand,  the  length 
should  be  chosen  only  large  enough  so  that  the  light  ray  at  the  end  of 
the  body  still  is  approximately  at  the  same  distance  from  the  wall  as 
at  the  initial  point  of  the  body. 
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3.  Test  Setup  and  Procedure"^" 

(a)  Plate-  -  The  plate  was  100  by  100  centimeters  and  1-0  centi- 
meter thick.   It  consisted  of  two  zinc  sheet  plates  1  millimeter  thick 
which  by  interposing  a  frame  of  8-millimeter  U-brass  were  soldered  onto 
each  other.   In  the  interspace  there  was  placed  a  heating  coil,  insu- 
lated by  asbestos  boards,  which  consisted  of  36  meters  of  chrome-nickel 
wire  of  1-millimeter  diameter  with  a  total  resistance  of  51  ohms.   In 
air  of  room  temperature,  the  plate  could  be  kept  at  100°  C  surface  tem- 
perature with  5.6  amperes.   The  plate  was  suspended  on  two  bicycle 
wheel  spokes  which  were  attached  to  the  upper  strip  of  the  frame. 
Through  the  latter  the  current  was  also  conducted  to  the  interior.   For 
measuring  the  surface  temperature  a  silver-constantan  thermocouple  was 
used  which  was  flatly  soldered  on  about  5  centimeters  from  the  upper 
rim  of  the  plate  and  the  wires  of  which  extended  a  few  more  centimeters 
quite  close  to  the  plate.   This  simple  temperature  measurement  could  be 
applied  because  extreme  accuracy  was  not  required.   Figure  20  shows 
the  plate  suspended  with  the  current -supplying  wires,  and  the  thermo- 
couple and  two  pendulums  each  with  two  spheres  which  served  to  determine 
the  vertical  and  made  possible  the  mutual  comparison  of  the  different 
photographs . 

(b)  Cylinder.  -  The  cylinder  consisted  of  two  layers  of  0.7- 
millimeter-thick  zinc  sheet  which  was  rolled  over  two  bicycle  wheel 
rims  and  soldered  to  their  upper  sheathing  strip.   Between  the  sheets 
there  was  placed  the  heating  coil  insulated  by  asbestos  board,  the  coil 
being  of  54  meters  of  chrome-nickel  wire  of  1-millimeter  diameter  with 
a  total  resistance  of  77  ohms,  which  brought  the  cylinder  in  air  of 
room  temperature  to  100°  C  surface  temperature  with  a  current  of  4.1 
amperes.   The  cylinder  was  suspended  on  two  bicycle  wheel  spokes.   The 
current  leads  and  thermocouple  were  similarly  located  on  the  upper 
sheathing  strip.   The  mean  length  of  the  cylinder  was  100.2  centimeters, 
and  its  mean  outer  diameter  58.45  centimeters  (fluctuating  between 
58.35  in  the  horizontal  and  58.75  over  the  soldered  seam).   The  frontal 
surfaces  were  closed  off  with  asbestos  board.   Figure  21  shows  a  view 
with  the  suspension,  current  leads,  and  the  two  pendulums. 

(c)  Setup.  -  Since  it  is  very  difficult  to  produce  a  parallel 
bundle  of  rays  of  the  diameter  of  the  body  here  employed  by  means  of 
lenses,  such  a  beam  was  replaced  by  the  light  from  a  strongly  screened 
(as  a  rule  of  3-mm  aperture)  arc  lamp  32  meters  from  the  heated  body. 
The  distance  between  the  body  and  the  screen  was  8  meters.   The  appa- 
ratus was  set  up  on  the  largest  floor,  of  42-meter  length,  of  the . 


^° Acknowledgment  is  made  to  Herr  Dr.  R.  Weise  and  Dr.  H.  Kurzweg, 
Leipzig.   For  the  test  setup,  acknowledgment  is  made  to  Master  Mechanic 
G.  Hentsch,  Leipzig. 
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Institute.   Because  of  the  window  in  one  of  the  walls  and  the  connec- 
tion to  the  staircase  halls,  it  was  not  possible  to  exclude  entirely 
currents  of  the  room  air  so  that  only  at  high  temperature  differences 
(above  50°  C)  could  useful  schlieren  photographs  be  obtained.   The 
existing  dissymmetry  of  the  critical  points  on  the  two  sides  of  the 
bodies,  both  for  the  plate  and  cylinder  (see  table  V)  was  due  to  a 
unilateral  convection  flow  from  the  window  wall.   The  horizontal 
adjustment  of  the  test  bodies  was  effected  by  means  of  a  water  balance; 
the  adjustment  parallel  to  the  light  cone,  by  measuring  the  shadow  of 
the  cold  body. 

(d)  Photography^'^.  -  The  schlieren  pictures  were  obtained  by  two 
methods.   In  one  method  the  shadow  pictures  were  taken  on  the  screen 
by  means  of  a  lens  and  camera  (13x18  cm,  Zeiss-Tessar  1:2.7).   In  spite 
of  highly  sensitive  plates  (Agf a-Superpan  16/10  DIN) ,  exposure  times 
of  5  to  25  seconds  were  required  on  account  of  the  low  brightness,  be- 
cause of  the  strong  screening  and  large  distance  from  the  arc  lamp. 
For  determining  a  time  mean  value  of  the  fluctuating  critical  points, 
these  time  pictures  were,  however,  particularly  well  suited.   The 
photograph  was  taken  either  somewhat  from  the  side  (plate)  or,  in  order 
to  avoid  distortions,  exactly  central  (partly  in  the  case  of  the  cyl- 
inder); the  camera  was  set  up  in  the  shadow  of  the  cylinder  so  that  the 
legs  of  the  tripod  threw  an  additional  shadow  (fig.  23).   The  second 
method  consisted  of  a  direct  illiimination  of  one  or  several  photoplates 
simultaneously  (Agfa-Isochron  13/l0  DIN,  10x15  cm),  which  were  brought 
to  the  position  of  the  screen.   The  greater  brightness  obtained  in  this 
way  made  possible  instantaneous  photographs  (through  light  flashes  from 
the  arc  lamp),  which  in  connection  with  their  greater  scale  provided 
a  view  in  instantaneous  detail  of  the  transition  from  laminar  to  turbu- 
lent flow. 


4.  Evaluation  of  the  Schlieren  Pictures 

Figures  22  and  23  show  time  exposures  of  the  plate  and  cylinder 
taken  by  means  of  lens  and  camera.   In  the  center  the  magnified  shadow 
of  the  body  by  the  thickness  of  the  temperature  boundary  layer,  the 
"inner  caustic  curve"  surrounding  it,  and  further  toward  the  outside 
the  "outer  caustic  curve"  arising  from  the  rays  near  the  wall.   The 
dissymmetry  of  the  lower  edge  of  the  plate  is  due  to  the  disturbances 
of  the  room  air  that  were  mentioned  previously.   The  upper  edge  of  the 
plate  is  immediately  below  the  frame  of  the  picture  in  the  "necking" 
of  the  central  shadow.   In  the  upper  stagnation  point  of  the  cylinder 
are  seen  the  shadows  of  the  suspension  wires,  current  leads,  and  thermo- 
couple wires;  below  are  seen  the  three  tripod  legs  for  the  central 
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photograph.   Figure  24  shows  for  simultaneously  ohtained  instantaneous 
photographs  of  the  left  °  side  of  the  plate  (in  the  case  of  direct 
illumination  ahout  25  cm  from  the  lower  edge  starting  with  small  inter- 
spaces of  about  1  cm),  a  region  of  about  50-centimeter  height:  at  24(a) 
and  the  lower  part  of  24(b),  laminar;  in  the  upper  part  of,  25(b),  a 
transition;  24(c)  and  24(d),  turbulent.   Figure  25  shows  three  simul- 
taneously taken  instantaneous  pictures  of  the  cylinder  with  direct 
illumination:  25(a)  at  the  lower  stagnation  point,  laminar;  25(b)  at 
the  left  azimuth  of  120°,  turbulent;  25(c)  at  the  upper  stagnation 
point,  warm  air  stream. 

As  shown  by  figures  22  to  25,  the  outer  caustic  curve  of  the  rays 
near  the  wall  is  particularly  suitable  for  determining  the  critical 
location.   If  this  line  is  broken  up,  then  turbulence  certainly  exists; 
whereas  the  flow  regions  farther  from  the  wall  are  more  easily  disrupted 
through  external  disturbances.   In  the  laminar  (lower)  region,  the 
caustic  curve  is  sharp  and  is  completely  stationary;  in  the  turbulent 
(upper)  region  it  is  torn  into  several  pieces  and  in  disordered  motion. 
In  the  time  photographs,  this  means  complete  washing  out.   The  critical 
position  itself  fluctuates  upwards  and  downwards,  at  one  time  very 
rapidly  and  at  another  slowly.   For  evaluating  the  time  photographs 
there  was  taken  as  the  "time  mean  value  of  the  critical  position"  the 
end  of  the  region  in  which  the  existence  of  an  outer  caustic  curve 
through  differences  in  brightness  could  still  be  determined  outwards 
and  inwards  (arrows).   Above  the  critical  position  the  heated  layer 
rises  as  a  whole  outwards  as  broad,  uniformly  bright  strips.   Comparison 
of  the  photographic  observation  thus  determined  with  the  mean  critical 
position  estimated  from  naked  eye  observation  over  a  period  of  time 
showed  good  agreement  and  therefore  the  justification  of  this  procedure. 

The  critical  positions  obtained  by  this  procedure  are  collected 
together  with  the  remaining  test  values  in  table  V.   Each  critical  point 
is  the  mean  value  of  four  determinations  two  of  which  were  made  with 
direct  illumination,  and  two  with  a  magnified  projection  picture.   The 
Gr  values  are  formed  in  the  case  of  the  cylinder  with  the  diameter,  in 
the  case  of  the  plate  with  the  critical  height  as  characteristic  length. 
By  means  of  equation  (74)  for  the  plate  and,  with  f(xj^j,),  equation  (84) 

for  the  cylinder,  the  critical  Reynolds  numbers  of  the  boundary  appear- 
ing in  the  last  column  were  computed.   Their  mean  values  of  303  for  the 
plate  and  285  for  the  cylinder  are  weighted  according  to  the  different 
illimination  times  of  the  individual  pictures.   From  the  mean  value 

Rs,   =  303  the  value  1.0x10  is  obtained  for  the  mean  critical  Gr 
kr 

for  the  plate. 


On  the  right  of  the  figure;  the  pictures  with  direct  ill\imination 
are  laterally  interchanged. 
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5.  Discussion  of  Results 

The  mean  value  of  303  for  the  critical  Reynolds  numher  of  the 

boundary  later  for  the  plate  and  285  for  the  cylinder  can  he  regarded 

in  agreement,  in  view  of  the  accuracy  of  determination  of  a  critical 

number  in  general  and  the  small  number  of  the  schlieren  pictures  and 

disturbances  of  the  room  air  in  particular.   This  is  understandable 

from  the  fact  that  the  velocity  profiles  of  both  flow  processes  at  any 

distance  from  the  lower  edge  or  lower  stagnation  point  are  affine  to 

each  other.   The  development  along  the  initial  stretch  occurs  in  a 

different  manner,  but  this  should  have  only  a  higher-order  effect.   As 

was  to  be  expected,  the  critical  number  of  the  free  profile  Re   *  300 

kr 

lies  considerably  below  the  critical  number  of  the  profile  in  the  case 
of  forced  flow  (plate),  which,  referred  to  the  displacement  thickness 
of  Burgers  and  Hansen,  amounts  to  950  in  the  wind  tunnel,  and  accord- 
ing to  more  recent  tests  in  the  Gottingen  water  tunnel  for  very  quiet 
approach,  it  amounts  to  about  1400.   It  would  be  of  interest  to  know 
whether  a  theoretical  stability  investigation  for  the  free  profile,  the 
analog  of  the  familiar  computations  in  the  forced-flow  case,  would  have 
to  be  conducted  in  order  to  give  this  increased  instability. 

With  the  determination  of  the  critical  number  for  the  vertical 
plate  and  the  horizontal  cylinder  and  reduction  to  the  critical  Reynolds 
number  of  the  boundary  layer,  the  data  of  other  investigators  will  be 
compared  with  regard  to  the  transition  from  laminar  to  turbulent  flow. 

E.  Schmidt  (ref.  26)  notes  that  according  to  tests  in  air,  the  flow 
at  about  100°  C  plate  temperature  remains  laminar  up  to  a  plate  depth 
of  about  50  centimeters,  corresponding  to  a  critical  Gr  of  8><10°,  in 
good  agreement  with  the  values  obtained  herein. 

M.  Jakob  and  W.  Linke  (ref.  19)  determined  from  a  combination  of 
the  measurements  of  various  investigators  on  vertical  plates,  vertical 
and  horizontal  cylinders,  block,  and  sphere,  as  the  transition  point 
from  the  GrV^  to  the  GrV^  law,  the  critical  number  of 
(Gr.Pr)j^^  3.0x10"^,  that  is,  for  diatomic  gases  Gr^^  «.  4.0x10"^.   The 

difference  as  compared  with  the  present  determinations  should  in  large 
part  be  ascribed  to  the  fact  that  hydrodynamically  unequivalent  bodies 
(Reynolds  number  of  the  boundary,  see  Part  III,  sec.  5a,  eq.  (88))  were 
here  treated  together.   If  account  is  taken  only  of  the  tests  on  the 
vertical  plates  and  vertical  cylinders  (Part  III,  sec.  6a),  which  for 
large  cylinder  diameters  and  large  Gr   (thin  boundary  layer)  are  at 
least  approximately  equivalent  hydrodynamically,  there  is  obtained  a 
critical  Gr  of  4x10^,  which  is  therefore  quite  close  to  the  present 
determinations . 
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E,  Griffith  and  A.  H.  Davis  (ref.  17)  investigated  the  local 
distrihution  of  the  heat  transfer  on  a  270-centimeter-high  wall,  which 
consisted  of  25  individual  elements.   They  obtained  a  minimum  value  of 
the  heat  transfer  at  ahout  38-centimeter  height  (their  own  determina- 
tion) which  they  regarded  as  the  transition  from  the  laminar  to  the 
turbulent  condition.   This  agrees  in  order  of  magnitude  with  the  pres- 
ent values  (columns  between  the  individual  elements),  where  it  should  be 
further  clarified  which  point  of  their  curve  is  to  be  regarded  as  the 
"critical  point"  proper. 


6.  Travel  of  the  Critical  Azimuth  for  the  Horizontal  Cylinder 

On  the  basis  of  the  given  theoretical  solution  of  the  flow  about 
the  horizontal  cylinder  and  the  point  of  transition  experimentally 
determined  for  a  Gr  ■=  10^,  the  travel  of  the  critical  azimuth  for  a 
variable  Gr  can  be  computed  if  it  is  assumed  that  the  turbulence 
always  occurs  if  the  boundary-layer  Re  has  the  same  value  of  285. 
This  reasonable  assumption  will  correspond  all  the  more  to  fact  since, 
as  has  just  been  shown,  even  for  the  vertical  plate  with  a  quite  differ- 
ent boundary-layer  development,  the  turbulence  is  set  up  at  the  same 
critical  number  of  the  boundary  layer. 

From  equation  (84),  the  following  relation  is  obtained,  with  the 

value  Re   =  285: 
°kr 

which  is  evaluated  in  table  VI  and  plotted  on  semilogarithmic  grid  in 
figure  26.   In  accordance  with  the  basic  considerations,  for  a  Gr  =  10^ 
the  start  of  turbulence  occurs  for  an  azimuth  of  about  120°.   Complete- 
ly laminar  flow  about  the  cylinder,  for  extension  of  the  theory  up  to 
the  upper  stagnation  point,  occurs  for  Gr<3^108.   The  transition  point 
reaches  the  equator  for  Gr  =  3X10  . 

For  greater  clarity,  the  Gr  will  be  expressed  for  cylinders  of 
varying  diameters,  which  are  in  air  of  20°  C  with  a  surface  temperature 
of  100°  C.   From  the  definition  of  Gr,  for  0  =  80  (°C), 
(3  =  1/293  (°C-^),  v(l00°)  =  0.231  (cm^/sec),  the  following  relation 
obtains : 

'Gr, 


d[ci§     =     f^xl03jl/3  (103) 


which  with  the  aid  of  the   critical     Gr,^     values   from  equation   (102)    has 

similarly  been  evaluated  in  table  VI  and  plotted  in  figure   27.      Corre- 
sponding to  the  basic   considerations,   the   critical  azimuth  for  the 
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relations  considered  for  a  cylinder  of  60-centimeter  diameter  is  120°. 
For  cylinders  below  41  centimeters,  no  turbulence  occiirs  at  all.   The 
transition  point  reaches  the  equator  for  a  cylinder  of  84-centimeter 
diameter. 

For  practical  applications,  the  result  is  thus  obtained  that  for 
heated  horizontal  conducting  pipes  up  to  40  centimeters  in  diameter  and 
surface  temperatures  of  100°  C  in  room-temperature  air  the  previously 
given  formulas  are  valid  for  the  heat  transfer  in  air  (in  particular, 
eqs.  (60)  and  (65)).   For  accurate  computation,  the  nondimensional 
formulation  of  figure  26  is  naturally  to  be  applied.   The  formulas  can, 
in  the  absence  of  other  suitable  data,  still  be  applied  up  diameters  of 
60  centimeters,  because  two  thirds  of  the  cylinder  periphery  is  still 
in  laminar  flow.   There  should  not  likewise  be  much  of  a  change  for 
considerably  higher  surface  temperatures  (200°  or  300°  C). 

The  experimental  determination  of  the  variation  of  the  critical 
azimuth  for  the  cylinder  with  the  Gr,  through  varying  the  temperature, 
remained  unsuccessful.   For  temperature  difference  up  to  about  50°  C 
the  outer  caustic  curve  is  still  not  sufficiently  far  from  the  central 
shadow  in  order  to  be  perceived  separately  from  the  latter,  as  is  re- 
quired for  determining  the  critical  number.   For  this,  the  distance 
between  the  cylinder  and  the  screen  would  have  to  be  considerably  in- 
creased, which  is  not  possible  on  account  of  space  requirements.   Going 
beyond  the  usually  employed  temperature  differences  of  80°  C,  because 
of  the  strong  increase  in  v,  no  longer  gives  an  increase  in  Gr  and 
moreover  gives  increasingly  stronger  deviations  from  the  assumption  of 
moderate  temperature  differences,  which  is  at  the  basis  of  the  theory. 
The  travel  of  the  critical  azimuth  with  Gr  is  best  determined  by 
varying  the  cylinder  diameter,  as  follows  from  figure  27. 


V.   SUMMARY 

From  the  numerous  tests  already  available  on  the  heat  transfer 
from  horizontal  pipes  and  wires  in  diatomic  gases,  the  dependence  of 
the  nondimensional  heat-transfer  number  Nu  on  the  Grashof  number  Gr 
and  the  temperature  coefficient  Te,  which  enters  as  a  further  non- 
dimensional  factor  at  large  temperature  differences,  is  determined.   The 
effect  of  Te  on  the  heat  transfer  is  quantitatively  determined  for  the 
first  time.   It  is  particularly  large  in  the  region  of  small  Gr  (lO~ 
to  10),  where  on  the  average  a  decrease  in  Nu  by  22  percent  for  in- 
crease of  Te  from  Te  =  0  to  Te  =  1  is  obtained.   The  experimentally 
found  dependence  of  Nu  on  Gr,  Pr,  and  Te  can  then  obtain  a  qualita- 
tive theoretical  explanation. 
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For  the  region  of  Gr  (about  10  to  3^10  )  in  which  the  heat  trans- 
fer is  limited  to  a  thin  (in  comparison  with  the  cylinder  diameter), 
heated  layer  with  laminar  flow,  the  velocity  and  temperature  fields  and 
the  heat  transfer  are  quantitatively  computed  from  the  boundary-layer 
differential  equations  without  any  additional  empirical  values;  these 
computations  are  found  in  good  agreement  with  the  available  measure- 
ments.  In  particular,  the  one-fourth  power  law  of  the  heat  transfer  is 
obtained  theoretically  as  Nu  =  0.37  GrV^. 

The  flow  and  heat-transfer  relations  thus  computed  (Re  of  the 
boundary-layer  flow,  variation  of  the  local  heat-transfer  coefficient 
along  the  cylinder  periphery,  mean  heat-transfer  coefficient,  and  total 
heat  given  off  by  the  cylinder)  are  compared  with  the  already  known 
relations  for  the  vertical  plate.   Among  other  results,  the  depth  of  a 
vertical  plate  is  determined  which  for  free  convection  shows  equal  flow 
condition  and  equal  total  heat  transfer  at  a  given  cylinder,  so  that  a 
simple  computation  is  made  possible  for  converting  the  heat-transfer 
data  for  a  plate  to  those  of  a  cylinder  and  conversely. 

In  order  to  know  the  upper  limit  of  validity  of  the  laminar -flow 
and  heat-transfer  computations  and  the  laminar  heat-transfer  formulas, 
the  start  of  turbulence  was  determined  by  schlieren  photographs  on  a 
vertical  plate  and  a  horizontal  cylinder  of  sufficient  size.   This 
occurs  for  different  values  of  Gr,  namely,  Gr  =  1. 0^*10^  for  the  plate 
and  Gr  =  3.5^10°  for  the  cylinder,  but  for  equal  Re  of  the  boundary- 
layer  flow,  namely,  for  Re  =  300.   In  this  way,  the  critical  Reynolds 
number  of  the  velocity  profile  of  the  free  flow  was  determined  for  the 
first  time;  as  compared  with  other  velocity  profiles,  this  profile  is 
marked  by  the  presence  of  a  maxim\im  value  and  a  point  of  inflection. 
Likewise  computed  with  the  Grashof  number  was  the  'critical  azimuth'  on 
the  cylinder  at  which  the  laminar  flow  passes  into  the  turbulent  flow. 


REFERENCES 

1.  Hermann,  E.:   Phys.  Z.,  Bd.  33,  1932,  p.  425. 

2.  Nusselt,  W.:  Gesundh.  Ing.,  Bd.  38,  1915,  p.  477. 

3.  Davis,  A.  H. :   Coll,  Res.  Nat.  Phys.  Lab.,  vol.  19,  1926,  p.  193. 

(See  also  Phil.  Mag,,  vol.  43,  1922,  p.  329.) 

4.  Grbber,  H.,  und  Erk,  S.  :   Die  Grundgesetze  der  Warmeiibertragung. 

Julius  Springer  Berlin,  1933. 

5.  Schmidt,  H. :   Ergebn.  exakt,  Naturwiss,,  Bd.  7,  1928,  p.  342. 


54  NACA  TM  1366 


6.  Grober,  H. :   WarmeulDertragung,  Julius  Springer  (Berlin),  1926. 

7.  Ayrton,  W.  E.,  and  Kilgour,  H. :   Phil.  Trans.  Roy.  Soc.  (London), 

vol,  A183,  1892,  p.  371. 

8.  Langmuir,  I.:   Phys.  Rev.,  vol.  34,  1912,  p.  401. 

9.  Langmuir,  I.:   Trans.  Am.  Inst.  Elec.  Eng.,  vol.  31,  1912,  p.  1229. 

(See  also  Proc.  Am.  Inst.  Elec.  Eng.,  vol.  32,  1913,  p.  407.) 

10.  V.  Bijlevelt,  J.  S.  :   Die  kiinstliche  Konvektion  an  elektrischen 

Hitzdrahte.  Diss.,  T.  H.  Dresden,  1915. 

11.  Kennelly,  A.  E.,  Wright,  C.  A.,  and  v.  Bijlevelt,  J.  S. :   Trans. 

Am.  Inst.  Elec.  Eng.,  vol.  28,  1909,  p.  363. 

12.  Petavel,  J.  E. :   Phil.  Trans.,  Bd.  (A)l97,  1901,  p.  229. 

13.  Wamsler,  F. :   VDI  Forschungsheft  98/99  (Berlin),  1911. 

14.  Nusselt,  W.  :   Diss.  T.  H.  Miinchen,  1907. 

15.  Koch,  W. :   Beih.  z.  Gesundh. -Ing.  Reihe  1,  Heft  22  (Berlin  und 

Munchen),  Oldenbourg  1927. 

16.  Jodltauer,  K. :   Forsch.  Ing.-Wes.,  Bd.  4,  1933,  p.  157. 

17.  Griffiths,  E.,  and  Davis,  A.  H. :   Transmission  of  heat  by  radiation 

and  convection.   (London),  1931. 

18.  King,  W.  J.:   Mech.  Eng.,  vol.  54,  1932,  p.  347. 

19.  Jakob,  M.,  and  Linke,  W. :   Forsch.  Ing.-Wes.,  Bd.  4,  1933,  p.  75. 

20.  Davis,  A.  H. :   Phil.  Mag.,  vol.  44,  1922,  p.  920. 

21.  Schmidt,  Ernst,  und  Beckmann,  Wilhelm:   Das  Temperatur-  und 

Geschwindigkeitfeld  vor  einer  Warme  abgebenden  senkrechter  Platte 
bei  naturlicher  Konvektion.   Tech.  Mech.  u.  Thermodynamik,  Bd.  1, 
Nr.  10,  Okt.  1930,  pp.  341-349;  cont.,  Bd.  1,  Nr.  11,  Nov.  1930, 
pp.  391-406. 

22.  Schmidt,  E. :   Forsch.  Ing.-Wes.,  Bd.  3,  1932,  p.  181. 

23.  Hermann,  R. :   Phys.  Z.,  Bd.  34,  1933,  p.  211. 


NACA  TM  1366  55 

24.  Beckmann^  W. :   Forsch.  Ing.-Wes.,  Bd.  Z,    1931^  p.  165. 

25.  Weise,  R. :   Forsch,  Ing.-Wes.,  Bd.  6,  1935,  p.  281. 

26.  Jakol),  M.,  und  Linke,  W.  :   Phys.  Z.,  Bd.  36,  1935,  p.  275. 

27.  Schmidt,  E. :   Z.  VDI,  Bd.  76,  1932,  p.  1025. 


Translated  by  S.  Reiss 
National  Advisory  Committee 
for  Aeronautics 


56 


NACA  TM  1366 


TABLE 

I.  -  SUMMARY  OF  TESTS  BY  DIFFERENT  INVESTIGATORS 

Investigator 

Ayrton, 
Kilgour 

Langmuir 

Bijlevelt 

Kennelly, 

Wright, 

Bylevelt 

Petavel 

Warns ler 

Koch 

Substance 

9  Pt  wires 

5  Pt  wires 

6  wires: 
Ta,  Pt,  Fe, 
Cu,  Ag,  Ni 

3  Cu  wires 

1  Pt  wire 

Pipes: 
6  wrought 
Fe;  6  Cu; 
1  cast  Fe 

4  steel 
pipes 

Cylinder  diam. , 
cm 

0.0031- 
0.0356 

0.00404- 
0.0510 

0.0043- 
0.1000 

0.01143- 
0.06907 

0.1106 

2.05-8.9 
5.9 
5.9 

1.4- 
10.05 

Temperature  of 
body,  t„,  °C 

40-300 

227-1027 

46-239 

— 

27.6-188.6 

Excess  tempera- 
ture, e,  °C 

15-180 

200-1000 

36-243 

13-174.4 

Room  tempera- 
ture, t^,  °C 

10.5-15.9 

27 

17-22 

20  (approxj 

16 

10-29 

14.2-22 

Medium 

Air 

Air 

Air 

Air 

Air,  Hg, 
°2 

Air 

Air 

Pressure, 

750 

750 

750 

120-1900 

0.12- 
160  atm 

715 

715.7-722.3 

Cylinder  length, 
cm 

32.5 

100 

35 

150  (approx) 

45 

300 

138-198 

Surrounding 
space 

Horizontal 
pipe: 
length, 
32.5  cm; 

diam. 
5.08  cm 

Base,  lOOx 
15  cm2; 
height, 
30  cm 

Height, 
300  cm 

Vertical 
tank: 
length, 
152  cm; 
height, 
660  cm 

Horizontal 
pipe: 
length, 
45  cm; 
diam. 
2.06  cm 

Base,  7.35 
sq  meters; 
height, 
210  cm 

Height, 
400  cm 

Ratio  of  space 

height  to 
cylinder  diam. 

1640: 143 

7430: 588 

70,000:3000 

58,200:9640 

18.6 

102:24 

286:40 

Te  number 

0.30-1.00 

0.67-3.33 

0.097-0.758 

0.052-0.61 

0.69-3.46 

0.13-0.82 

0.045- 
0.606 

Relative 
decrease  of 
Nu  for 
ATe  =  1, 
percent 

7.6-16.8 

7.2-2.2 

8.4-45.6 

25.2 

Irregular 

None 

None 

Mean  value 
of  Ku  for 
ATe  =  1, 
percent 

13.2 

Systematic 

25.6 

25.2 





Figure 
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TABLE  V.  -  EVALUATION  OF  SCHLIEREN  PHOTOGRAPHS  FOR  OCCURRENCE 

OF  TURBULENCE 
(a)  Vertical  plate. 


Number 

mm 

too, 

°C 

Exposure 
time, 
sec 

°C 

Critical  height,   h      ,    cm 

^^hkr 

^6kr 

Left 

Right 

Middle 

24 
25 
26 
27 

748 
748 
748 
748 

16.5 
16.5 
16.5 
16.5 

20 
20 
25 
20 

78 
113 
100 
100 

76.0 
53.8 
66.2 
60.9 

55.9 
47.3 
60.8 
47.3 

66.0 
50.6 
63.5 
54.1 

13.42  xlO^ 
6.83  xlO^ 

13.13  xlO^ 
8.12  xIqQ 

324 
274 
323 
286 

Mean:   303 


(b)  Horizontal  cylinder.   Diameter,  d,  58.45  centimeters, 


Number 

b, 
mm 

°C 

Exposure 
time, 

tv^ 
o^ 

o 

Critical  azimuth,    x. 

G^hkr 

'^  kr' 

^6kr 

Left 

Right 

Middle 

sec 

6 

760 

19 

10 

102 

125.5° 

o 
123.5 

124.5° 

10.26x10^ 

4.55 

291 

10 

760 

19 

20 

102 

125.8° 

110.0° 

117.9° 

10.26  xIQS 

4.35 

278 

11 

760 

19 

5 

102 

120.9° 

116.0° 

118.5° 

10.26x10^ 

4.37 

279 

12 

760 

19 

20 

102 

126.7° 

117.6° 

122.1° 

10.26  xl08 

4,48 

286 

17 

752 

18 

20 

102 

123.8° 

115.0° 

119.4° 

10.23  xlO^ 

4.39 

281 

18 

752 

18 

5 

102 

144.5° 

129.5° 

137.0° 

10.23  xIqQ 

4.90 

313 

Mean:   285 


TABLE  VI.  -  TRAVEL  OF  TURBULENCE  TRANSITION  POINT  (CRITICAL  AZIMUTH  xj^^,) 
©ER  OF  CYLINDER   Gr^  OR  WITH  CYLINI 
[t^,  100°  C;  toe,  20°  C;  b,  760  mm.] 


WITH  GRASHOF  NUMBER  OF  CYLINDER   Gr^  OR  WITH  CYLINDER  DIAMETER   d 


^kr 

180° 

165° 

150° 

120° 

90° 

60° 

30° 

15° 

10° 

5° 

G^d 
d,    cm 

3,52 
xIqS 
41.3 

4.30 
xIqS 

44.1 

5.42 
xIqS 

47.7 

1.07 
xl09 
59.9 

2.96 
xl09 
84.0 

1.37 

xlolO 

140 

2.10 

xlO^l 

347 

3.35 
xl0l2 

890 

1.69 

xl0l2 

1500 

2.72 

xlOl^ 

3790 
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log  Gr 
Figure  1.  -  Convection  tests  of  Ayrton  and  Kllgour  on  nine  platinum  wires. 
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Figure  2.  -  Convection  tests  of  Langmuir  on  five  platinum  wires. 
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Figure  3.  -  Convection  tests  of  Bljlevelt  on  six  wires  of  tantalum,  platinum 

iron,  copper,  silver,  and  nickel. 
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Figure  4.  -  Convection  tests  of  Kennelly,  Wright,  and  Bijlevelt  on  three 
copper  wires. 
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Figure  5.  -  Convection  tests  of  Petavel  on  one  platinum  wire. 
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Figure  6.  -  Convection  tests  of  Wamsler  on  one  copper,  one  cast  iron,  and  six 
wrought  iron  pipes. 


62 


MCA  TM  1366 


log  Gr 
Figure  7.  -  Convection  tests  of  Koch  on  four  steel  pipes. 
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(extrapolated).   Limiting  law  of  small  temperature  differences. 
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Figure  9.  -  Comparison  of  test  results  of  various  Investigators 
for  Te  =  0.65. 
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Figure  10.  -  Summary  of  all  test  results  with  diatomic  gases  (Pr  =  0.74)  with 
the  adjustment  curve  giving  the  most  probable  heat-transfer  law  based  on 
present  information. 
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Figure  13.  -  Velocity  profiles  near  cylinder  referred  to  nondimensional  coordinates. 
Continuous  curve  represents  theoretical  solution.  Points  according  to 


measurements  of  Jodlbauer  for  2r  =  9  centimeters;  t^ 
Gr  =  3.76xl05. 
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Figure  14.  -  Temperature  profiles  near  cylinder  referred  to  nondimensional 
coordinates.  Continuous  curve  represents  theoretical  solution.  Points 
according  to  measurements  of  Jodlbauer  for  2r  =  9  centimeters;  ty  =  99.2° 
t^=  18.1°  C;  Gr  =  S.TGxIqS. 
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Figure  15.  -  Velocity  profiles  near  cylinder.  Continuous  curve  gives 

theoretical  solution.  Points  according  to  measurements  of  Jodlbauer  for 
2r  =  5  centimeters;  t„  =  104.6°  C;  tco=  18.1°  C;  Gr  =  6.54xlo5. 
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Figure  16.  -  Temperature  profiles  near  cylinder.  Continuous  curve  gives 
theoretical  solution.  Points  according  to  measurements  of  Jodlbauer  for 
2r  =  5  centimeters;  t„  =  104.6°  C;  t„  =  18.1°  C;  Gr  =  6.54x10-. 
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Figure  17 .  -  Azimuth  function  g(x)  giving  variation  of  local  heat-transfer 

coefficient  or  of  a  reciprocal  distance  from  the  wall  which  is  characteristic 
for  profile  (for  example,  boundary- layer  thickness)  along  cylinder  perimeter. 
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Figure  18.  -  Comparison  of  heat  transfer  for  rectangular  plate  and  horizontal 
cylinder.  For  H  =  d,  under  otherwise  equal  conditions,  the  two  diagrams 
give  dimensional  heat-transfer  coefficients  a    to  scale. 
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Figure  19.  -  Ratio  of  total  heat  transfer  of  horizontal  cylinder  and  vertical 
plates  (on  both  sides)  of  different  heights  under  otherwise  equal  conditions. 
Ratio  of  boundary- layer  Re  at  upper  stagnation  point  of  cylinder  and  upper 
edge  of  plate.  Abscissa,  heat  transfer  and  Re;  ordinate,  plate  height. 


Figure  20.  -  Vertical  plate  100x100x1  centimeters  for  determining  occurrence 
of  turbulence  for  free  convection. 
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Figure  21.  -  Horizontal  cylinder  (58.5-cm  diani.j 
occurrence  of  turbulence  for  free  convection. 


100- cm  length)  for  determining 


Figure  22.  -  Schlleren  photograph  of  heat  transfer  at  the  vertical  plate, 
photograph  (25  sec)  with  lens  and  camera.  Surface  temperature,  t__ 
t„=  16.5°  C;  b  =  748  millimeters;  critical  height  (arrows),  left"66.2 
centimeters,  right  60  centimeters;  critical  Grashof  number,  13.1x10  ; 
critical  boundary-layer  Ee, 323. 
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i^u.z'e  23 .  -  SciLliere-  piotograph.  of  tie  Lea  -  "ransf  er  ar  ~hB  h.c 
Time  photographs  ("0  sec)  vith  lens  and  camera,  ^t-,-,  102°  C;  1 
b,  752  atilliaieters;  (ir  of  the  cylinder,  ICSiLC^;  critical  a: 
left  125.8°,  ri^t  115. 0°;  critical  caimdarj- layer  3e,   2£1. 
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r-LTS  24.  -  Schllerei:  pLot-ogrspIis  of  best:  transfer  at  one  side  of  vertical 
;l£te.     5l3iultareoTi£  Ir^tantareo-us  pbstogrsphs     (a'cc«t  1/2C  sec)   vlth  Ir^dlrect 


— ---•   cio-^  -; 


=<>  - 


part  of  (b),  traiisitiOE.;   (c)   and  (d),  turbulent. 
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Figure  25.  -  Schlieren  photographs  of  heat  transfer  at  the  horizontal  cylinder. 
Simultaneous  instantaneous  photographs  (l/20  sec)  with  Indirect  illumination, 
t^  =  102°  C;  toD=  18°  C;  b  =  752  millimeters;  Gr,  10.2x10^.   (a)  lower  stag- 
nation point,  laminar;  (b)  azimuth  of  about  120°,  turbulent;  (c)  upper  stagnation 
point,  upflow  of  warm  air  with  current  lead  wires  and  so  forth. 
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Figure  26.  -  Variation  of  position  of  start  of  turbulence  with  Gr 
number  of  cylinder.  Computed  for  critical  boundary- layer  Ee,  285. 
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Figure  27-  -  Variation  of  position  of  start  of  turbulence  Vith  the  cylinder 

diameter  for  surface  temperature  of  100°  C  in  air  at  20°  C  and  760  millimeters 
of  mercury.  Computed  for  critical  boundary- layer  Ee,  285. 
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